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INTRODUCTION: 


Year  Four  of  this  applied  research  program  continues  to  build  on  advances  by  our  collaboration 
in  the  development  of  novel  methods  and  advanced  hardware  for  low-magnetic-field  MRI. 
Without  major  innovation,  high-field  MRI  instruments  offer  limited  utility  for  imaging  TBI  in 
widely  deployable  contexts.  We  focus  our  research  effort  on  the  high-risk  and  critical  challenges 
that  must  be  solved  to  enable  deployment  of  a  transportable  human-head  MRI  system  applicable 
to  TBI  imaging  in  battlefield  medical  facilities.  Our  goal  is  to  establish  proof-of-principle  of  a 
suite  of  techniques  and  technologies  to  advise  future  development  of  a  field- deployable  device 
with  high  diagnostic  impact.  This  research  effort  has  two  specific  aims: 

Specific  Aim  1:  Develop  a  low- field  human-head  MRI  system  (LFI)  suitable  for  high- resolution 
multi-nuclear  imaging,  and  improve  the  ability  to  attain  brain  images  based  on  the  intrinsic  in 
vivo  'H  NMR  signal  in  this  scanner. 

This  includes  the  development  of  two  robust  low-field  scanner  hardware  methodologies  (elec¬ 
tromagnet  based  and  permanent  magnet  based),  the  development  of  novel  high-speed  parallel 
imaging  detection  systems,  and  work  on  advanced  adaptive  reconstruction  methods  including 
navigators  and  sparse  sampling. 

Specific  Aim  2:  Develop  injury- sensitive  MRI  based  on  converting  the  electron  spin  of  free  radi¬ 
cals  associated  with  injury  (specifically  TBI)  into  nuclear  polarization  using  the  Overhauser  effect 
and  subsequently  imaging  that  modified  nuclear  polarization  using  low- field  MRI  (OMRI).  Suc¬ 
cessful  demonstration  of  OMRI  of  free  radicals  associated  with  injury  will  be  directly  applicable 
to  the  MRI  systems  of  Aim  1,  enhancing  image-based  injury  specificity  and/or  shortening  scan 
acquisition  time. 

The  development  of  this  new  MRI  contrast  mechanism  may  provide  an  unambiguous  non- 
invasive  in  vivo  maker  for  cerebral  injury,  and  has  potential  for  assisting  the  imaging  of  TBI  at 
both  low  and  high  magnetic  fields. 
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BODY: 


Image  quality  in  the  6.5  mT  electromagnet  LFI  (Aim  1)  continues  to  improve  as  our  suite  of 
hardware  and  image  acquisition  strategies  matures.  Validation  of  image  quality  in  head-shaped 
phantoms  continues,  including  images  acquired  using  a  new  single-channel  receiver  coil  opti¬ 
mized  for  human  heads.  This  receiver  coil  design  demonstrates  performance  exceeding  that  of 
the  more  complex  multi-channel  designs  described  previously.  We  report  details  of  this  design, 
as  well  as  images  acquired  in  a  new  multi- compartment  head- shaped  phantom  developed  for  this 
work. 

Additionally,  we  continue  to  make  significant  progress  in  our  permanent  magnet  Halbach  LFI 
experiment  (Aim  1)  toward  the  acquisition  of  true  3D  datasets  in  the  Halbach  magnet  portable 
imaging  prototype  [1].  For  the  first  time,  three  key  ingredients  -  “built-in”  2D  encoding  using  a 
rotating  B0  encoding  field,  phase  encoding  with  a  linear  phase  RF  transmit  coil,  and  broadband 
RF  pulses  -  were  brought  together  to  enable  three-axis  encoding  in  an  inhomogeneous  B0  field 
with  no  gradient  coils. 

Our  previous  work  on  the  detection  of  free  radicals  (Aim  2)  focused  on  decreasing  the  image  ac¬ 
quisition  time  in  Overhauser  MRI  and  on  improving  the  detection  limit  to  low  concentrations  of 
free  radicals  in  synthetic  radical  solutions.  Work  on  this  task  continues  in  collaboration  with  Pro¬ 
fessor  of  Chemistry  D.  Whitney  King  (Colby  College.  Dr.  King’s  laboratory  has  developed  a  port¬ 
able  calibrated  superoxide  generator,  and  will  allow  us  to  directly  investigate  OMRI  sensitivity  to 
biologically  relevant  superoxide  at  nanomolar  to  micromolar  concentrations.  This  is  a  critical 
step  to  understand  the  biological  relevance  and  applicability  of  the  new  method  of  free  radical 
imaging  developed  under  this  grant.  We  determined  that  significant  improvements  were  still 
needed  in  the  manner  in  which  our  OMRI  detection  is  performed  before  we  could  attain  sensitiv¬ 
ity  to  superoxide,  and  to  that  end  have  improved  our  OMRI  detector  hardware  and  attained  a  3- 
fold  improvement  in  SNR  and,  critically,  wide-range  tunability  has  been  attained. 

Low  Field  Imagers  (LFIs) 

TASK  1A:  Low-field  MRI  Hardware  Development 

We  continue  to  optimize  hardware  development  for  the  human  head  LFI  test  bed  systems  (Aim 
1)  and  for  the  OMRI  system  (Aim  2). 
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The  electromagnet  LFI  (shown  below)  is  an  ideal  state- 
of-the  art  test  bed  for  all  of  the  novel  acquisition,  detec¬ 
tion  methodologies,  and  reconstruction  algorithms  in¬ 
cluding  navigators  and  sparse  sampling,  and 
additionally  will  provide  necessary  experience  and  data 
to  advise  optimal  construction  and  magnetic  field  for 
any  future  electromagnet-based  deployable  systems. 
Currently,  this  state-of-the-art  scanner  enables  high- 
performance  spectroscopy  and  8-channel  imaging  at  6.5 
mT,  and  is  fully  equipped  for  Overhauser  DNP  experi¬ 
ments. 

The  permanent  magnet  system  is  a  lightweight  (45  kg) 
and  portable  Halbach  array.  This  Halbach  array  scanner 
is  a  highly  specialized  and  potentially  disruptive  tech¬ 
nology  scanner  that  could  greatly  ease  both  the  cost  and 
burden  of  a  field-forward  instrument  purpose-built  for 
TBI  imaging.  This  Halbach  imager  contrasts  markedly 
with  the  electromagnet  LFI  in  that  it  has  a  highly  inho¬ 
mogeneous  magnetic  field,  but  we  use  this  inhomogene¬ 
ity  to  our  advantage  as  an  encoding  field  to  acquire  head  images  without  the  use  of  an  additional 
gradient  set. 

Halbach  array  (permanent  magnet)  LFI 

This  test  bed  scanner  is  based  around  a  lightweight  array  of  permanent  magnets  in  a  so-called 
Halbach  configuration.  This  magnet  is  ideal  for  portable  MRI  in  that  it  creates  a  relatively  uni¬ 
form  field  transverse  to  the  head  without  the  use  of  a  cryostat  or  power  supplies.  This  Halbach 
array  scanner  is  a  highly  specialized  scanner  and  a  potentially  disruptive  technology  that  could 
greatly  ease  both  the  cost  and  burden  of  a  field-forward  instrument  purpose-built  for  TBI  imag¬ 
ing.  This  Halbach  imager  contrasts  markedly  with  the  electromagnetic  LFI  in  that  it  has  a  highly 
inhomogeneous  magnetic  field,  but  we  use  this  inhomogeneity  to  our  advantage  and  use  it  to  ac¬ 
quire  head  images  without  the  use  of  an  additional  gradient  set. 

A  truly  portable  MR  system  has  the  potential  to  quickly  detect  brain  injury  at  the  site  of  injury. 
For  example  hemorrhage  detection  is  critical  for  both  stroke  patients  and  traumatic  brain  injury 
victims.  In  stroke,  rapid  distinction  between  a  hemorrhagic  and  non-hemorrhagic  event  could 
allow  administration  of  a  clot-busting  drug  such  as  tPA  (tissue  plasminogen  activator)  in  an  am¬ 
bulance  prior  to  transportation  to  the  hospital,  perhaps  advancing  this  time-sensitive  treatment 
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by  up  to  an  hour.  Subdural  hemorrhage  (or  hematoma)  is  a  form  of  traumatic  brain  injury,  in 
which  blood  gathers  between  the  dura  and  arachnoid  mater  (in  meningeal  layer)  and  is  likely  to 
be  visualized  on  course  resolution  (e.g.  5  mm)  T1  images. 


In  Y1  we  designed  and  built  the  very  portable  45  kg  Halbach  array  magnet  for  portable  MRI,  and 
mapped  the  magnetic  field  and  drift  using  NMR  measurements  at  3.3  MHz.  In  Y2  we  focused  on 
the  hardware  and  encoding  techniques  needed  to  turn  this  magnet  into  an  imager.  In  Y3,  we 
solved  some  critical  challenges  including  the  ability  to  acquire  2D  multi-channel  images.  Re¬ 
search  in  2014  (Y4)  saw  substantial  progress  toward  the  acquisition  of  true  3D  datasets  in  the 
Halbach  magnet  portable  imaging  prototype  [1].  For  the  first  time,  three  key  ingredients  -  “built- 
in”  2D  encoding  using  a  rotating  B0  encoding  field,  phase  encoding  with  a  linear  phase  RF 
transmit  coil,  and  broadband  RF  pulses  -  were  brought  together  to  enable  three-axis  encoding  in 
an  inhomogeneous  B0  field  with  no  gradient  coils. 

Halbach  array  LFI:  extension  to  3D  encoding 

Previously  in  2013,  we  achieved  2D  axial-plane  imaging  without  gradient  coils  by  using  a  rotat¬ 
ing,  45  kg,  77  mT  inhomogeneous  Halbach  cylinder  based  magnet.  The  built-in  field  variation  of 
the  permanent  magnet  array  was  used  as  a  rotating  spatial  encoding  magnetic  field  (rSEM).  As 
the  magnet  rotates  around  the  sample,  generalized  projections  onto  the  non-linear  SEM  are  ac¬ 
quired  as  a  spin  echo  train.  Images  (Figure  1)  are  reconstructed  from  the  projection  data  by  ap¬ 
plying  an  iterative  solver  to  the  full  encoding  matrix  calculated  using  the  signal  equation  [1], 


*6*  *8* 


Figure  1:  Ground  truth  phantom  with  text  (left)  contained  in  a  1.5cm  thick  cylinder  slab.  2D  im¬ 
aging  based  on  Halbach  rSEM  data  is  shown  at  right.  Parallel  receive  coils  are  used  to  disambigu¬ 
ate  the  non-bijective  mapping  of  the  nonlinear  encoding  field. 
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TRASE  encoding 


To  enable  imaging  in  the  3rd  direction  (along  the  axis  of  the  cylindrical  magnet,  see  Figure  2),  we 
use  the  Transmit  Array  Spatial  Encoding  (TRASE)  method  [2],  [3]  TRASE  is  a  BE  encoding 
method  that  requires  2  different  switchable  BE  phase  gradients  (often  equal  and  opposite)  along 
the  encoding  direction.  The  ID  TRASE  sequence  (Figure  3)  is  a  modified  RARE  spin-echo  train 
in  which  the  slope  of  the  Bf  phase  gradient  switches  for  consecutive  refocusing  pulses  (indicated 
by  the  +  or  -).  The  two  phase  gradients  correspond  to  two  k- space  origins  {-ki  and  ki),  and  the 
previous  k- space  point  is  flipped  across  the  k- space  origin  of  the  refocusing  pulse  to  traverse  k- 
space.  With  each  successive  pulse  in  the  echo  train,  the  object  is  modulated  by  additional  Aki  un¬ 
til  sufficient  spatial  encoding  is  achieved. 


#biW  =  2nk1x 


Figure  2:  (a.)  Y-Z  encoding  is  performed  with  the  rotating  SEM  method,  (b.)  TRASE  is  imple¬ 
mented  in  the  X  direction  using  a  switchable  linear  B1+  phase  along  X. 

Coil  hardware 

A  TRASE  array  was  designed  that  consists  of  two  nested  cylindrical  coils  (Figure  4).  Coil  1  is  a 
short  4-turn  birdcage  coil  (12  rungs,  18cm  diam.,  22cm  length)  [4]  that  produces  a  BE  field  in  Y 
with  a  cosine  shape  along  X,  B+y  (x)  =  |5+xy  |cos  (2n  k1x).  Coil  2  is  a  10  turn  Maxwell  coil 
(22cm  diam.,  18cm  length)  that  produces  a  BE  field  in  X  with  a  sine  shape  along  X,  B+X(x')  = 
±|£ixy|sin  (2n  ktx). 


8 


TRASE  n/J  X  rn 

RF  rH  \ 

n 

h  7t 

SEM 

Acquisition 

m 

warn 

-ki  kx 

k- space - 1 - • - 

0 

Figure  3:  The  TRASE  sequence  is  repeated  for  every  magnet  rotation  of  the  imaging  sequence. 

The  +/-  indicates  the  sign  of  the  TRASE  phase  slope. 

The  coils  are  tuned  to  the  Larmor  frequency  (3.29  MHz)  and  decoupled  with  a  transformer  (-20 
dB  isolation).  Figure  4b  shows  the  relative  magnitude  of  Blx  and  B,y  from  the  2  coils  measured 
with  a  pick-up  loop.  When  the  two  coils  are  driven  with  equal  \B+xy  |  magnitude  and  spatial  fre¬ 
quency,  kb  they  produce  the  desired  Bt+  field  variation  along  X,  By  (x)  =  \B+xy \e+l  2nk^x  (uni¬ 
form  magnitude  and  linear  phase).  A  switchable  180°  phase  shifter  is  added  in  the  Maxwell  coil 
RF  path,  which  applies  a  sign  change  to  every  other  refocusing  pulse.  The  pulses  with  the  phase 


Figure  4:  The  TRASE  array  is  a  nested  birdcage  and  Maxwell  coil,  b:  The  measured  magnitude 
from  Blx  of  the  Maxwell  coil  and  B  j}.  of  the  birdcage  coil  is  plotted 

Acquisition 

The  2D  rSEM  method  and  TRASE  are  performed  simultaneously  for  3D  imaging.  The  TRASE 
array  must  rotate  with  the  magnet  because  the  Y  direction  changes  as  the  magnet  rotates.  For 
every  magnet  rotation,  the  projection  onto  the  SEM  is  acquired  as  a  TRASE-modulated  16-echo 
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spin-echo  train  (each  echo:  readout  =  196  pt,  BW  =  40  KHz,  echo  spacing  =  9.7  ms).  Data  is  ac¬ 
quired  at  180  rotations  of  the  magnet  spaced  1°  apart.  Frequency-swept  WURST  pulses  are  used 
to  achieve  the  same  flip  angle  across  the  inhomogeneous  field  (BW  =  25  KHz,  [pw90,  pwl80]  = 
[6,3]  ms)  [5].  The  birdcage  coil  of  the  TRASE  array  is  used  as  the  receiver  coil. 

Correction  of  quadratic  phase  from  frequency-swept  pulses 

It  is  difficult  to  perform  TRASE  encoding  using  short,  rectangular  “hard”  pulses  in  the  presence 
of  the  Halbach  inhomogeneous  AB0  SEM  due  to  practical  limitations  on  the  available  RF  power, 
and  thus  the  minimum  achievable  refocusing  pulse  duration.  The  TRASE  point  spread  function 
degrades  rapidly  for  refocusing  pulse  angles  outside  the  range  150-195°  due  to  the  accumulation 
of  flip  angle  and  phase  errors  along  the  echo  train  [6] .  For  the  hard  pulses  used  in  early  published 
TRASE  images  [2]  the  refocusing  angle  falls  below  150°  when  AB0  >  40%  of  the  pulse  Rabi  fre¬ 
quency. 

To  overcome  this  limitation,  we  perform  TRASE  in  the  presence  of  strong  off-resonance  using 
WURST-40  pulses  [7]  with  phase  cycling  [5]  for  accurate  spin  excitation  and  refocusing  across 
also  isochromats.  Pulses  with  a  linear  frequency  sweep  impart  a  quadratic  phase  across  the  spin 
bandwidth,  but  if  the  pulse  sweep  rate  satisfies  the  condition  Rexc  =  Vi  Rref  and  the  AB0  SEM  is  con¬ 
stant  over  time,  then  the  refocusing  pulse  removes  the  phase  applied  by  the  excitation  pulse  and 
all  isochromats  refocus  the  echo  time  [8].  Further  readouts  in  the  echo  train  alternate  between 
echoes  and  “spectral”  echoes  [5]  that  resemble  the  frequency  distribution  of  spins  in  the  AB0 
SEM.  This  additional  RF  phase  modulation  must  be  removed  to  make  the  spectral  echoes  com¬ 
patible  with  TRASE. 

The  spectral  echo  signal  can  be  modeled  as  a  convolution  of  the  conventional  spin  echo  signal 
with  a  chirp  function,  exp(i27r/?t2),where  f  depends  on  the  sweep  rate.  Deconvolution  is  per¬ 
formed  by  either  (a.)  convolving  the  spectral  echo  with  a  scaled,  conjugated  version  of  the  chirp 
function,  or  equivalently  (b.)  in  the  Fourier  transform  domain  by  dividing  the  FFT  of  the  echo 
(ID  projection)  by  the  FFT  of  the  chirp  function  [9].  Since  the  FFT’ed  chirp  kernel  is  a  pure 
phase  term,  method  (b.)  requires  simply  taking  the  difference  between  the  phase  of  the  first  and 
second  projections  from  the  echo  train  (Figure  5).  This  phase  difference  is  used  to  correct  all 
subsequent  spectral  echoes. 
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Figure  5:  WURST  pulse  sequence  diagram  with  four  echoes  acquired  on  a  water  bottle  phantom. 
The  bracketed  block  is  repeated  as  needed.  Conventional  echoes  alternate  with  “spectral”  echoes 
modulated  by  the  quadratic  phase  of  the  frequency  swept  RF  pulses.  This  effect  is  compensated  by 
the  removal  of  the  quadratic  phase  from  the  projections  corresponding  to  each  spectral  echo.  Af¬ 
ter  phase  correction,  the  echo  real  and  imaginary  parts  overlay  for  all  echoes  in  the  train.  Phase 
cycling  scheme  is  as  follows:  <pexc  =  [90,  270,  180,  0]°,  (pref1  =  [0,  0,  90,  90 ]°,  0re/2  = 
[90,  90,  0,  0]°,  (pacq  =  [0,  180,  90,  270]°.  The  3rd  and  4th  echoes  in  the  block  are  multiplied 
by  tc  to  compensate  for  phase  cycling. 


To  validate  the  phase  correction  method,  a  ID  multi-tube  phantom  is  placed  in  the  Halbach  AB0 
SEM  and  a  train  of  echoes  is  acquired  using  TRASE.  WURST-40  excitation  and  refocusing  pulse 
durations  are  [ Texc ,  Tref]  =  [6,  3]  ms  and  sweep  rates  are  [Rexc,  Rref ]  =  [4.15,  8.3]  MHz/s  over  a 
25  KHz  bandwidth. 


Fourier  domain  phase  subtraction  faithfully  restores  spectral  echoes  so  that  their  real  and  imagi¬ 
nary  components  align  closely  with  those  of  conventional  echoes  (Figure  6).  The  TRASE  acquisi¬ 
tion  successfully  generates  projections  of  the  phantom  even  though  its  bandwidth  exceeds  10 
KHz  in  the  AB0  SEM.  By  contrast,  recognizable  projections  could  not  be  obtained  using  hard 
pulses  since  available  RF  power  constrained  the  minimum  pulse  duration  to  [75,150]  [ is  for 
[90°,  180°]  pulses,  providing  inadequate  bandwidth. 
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TRASE  k-space  1-D  projection 
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Figure  6:  The  first  8  corrected  echoes  (a.)  from  an  echo  train  shown  with  (red)  and  without 
TRASE  phase  modulation.  Orange  and  green  dots  at  each  echo  time  denote  positive  and  negative 
k-space  points,  respectively,  which  are  reordered  into  ID  k-space  (b.)  The  Fourier  transform 
yields  a  ID  projection  (c.)  of  the  water-filled  tube  phantom  (d.). 


3D  Reconstruction 


Once  the  raw  data  are  phase  corrected  as  described  above,  image  reconstruction  proceeds.  Step 
one  separates  the  echo  train  data  along  X  to  obtain  16  different  YZ-datasets  corresponding  to 
each  slice.  For  each  point  in  the  readout  window,  the  set  of  16  points  along  the  echo  train  popu¬ 
lates  a  k-space  vector  (Figure  3).  The  points  are  rearranged  and  FFT’ed  to  generate  a  ID  projec¬ 
tion  along  X  (Figure  6  and  Figure  7).  Only  the  k-space  lines  formed  from  the  center  points  of  the 
echoes  result  in  ID  projections  that  are  not  modulated  by  the  B0  SEM.  This  is  repeated  for  the 
echo  trains  from  each  of  the  magnet  rotations.  Step  two  reconstructs  each  set  of  echoes  (YZ- 
data)  into  a  2D  image  of  the  corresponding  slice.  The  encoding  matrix  for  each  slice  is  calculated 
using  2D  B0  field  maps  and  then  iterative  matrix  inversion  is  performed  to  separately  obtain  each 
2D  image. 

Results: 


Figure  7  shows  experimental  3D  imaging  results  using  a  7.5  cm  thick  phantom  with  1.5  cm  wa¬ 
ter-filled  compartments  spaced  3cm  apart  in  X.  A  ID  projection  along  X  shows  the  water-filled 
compartments.  Three  slices  are  reconstructed  from  data  at  X  =  [-3cm,  Ocm,  3cm]  using  the  cor¬ 
responding  2D  AB0  field  maps. 
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Discussion 


Figure  7:  Image  of  a  3D  phantom  consisting  of  three  water-filled  compartments  spaced  3  cm  apart 
in  X.  A  ID  projection  along  X  for  one  rotation  angle  is  shown.  Three  image  slices  were  recon¬ 
structed  using  the  appropriate  2D  field  maps. 


The  non-bijectivity  of  the  SEMs  results  in  image  aliasing  across  the  center.  This  is  usually  re¬ 
solved  using  a  surface  coil  array  [10],  but  as  a  first  step,  the  birdcage  coil  of  the  array  is  used  for  a 
receiver  coil  and  only  the  reduced  FOV  image  is  shown.  The  blurring  in  the  image  slice  is  due  to 
a  combination  of  the  spatially  varying  resolution  of  the  SEMs  and  systematic  calibration  errors  in 
our  model-based  reconstruction. 


We  have  demonstrated  proof-of-concept  3D  images  in  the  inhomogeneous  45  kg  Halbach  mag¬ 
net  without  the  use  of  gradient  coils.  TRASE  encoding  along  X  is  ideal  because  like  our  2D  meth¬ 
od  it  uses  spin  echo  trains,  and  it  avoids  the  need  for  a  B0  gradient  coil  and  gradient  power 
amplifier,  which  would  add  weight,  cost,  and  power  requirements  to  the  portable,  low-cost  sys¬ 
tem. 

We  have  further  shown  that  robust  TRASE  spatial  encoding  in  an  inhomogeneous  B0  field  can  be 
achieved  using  frequency-swept  WURST  pulse  echo  trains  when  a  simple  correction  method  is 
used  to  remove  the  unwanted  quadratic  phase  from  alternate  echoes. 
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High-efficiency  helmet  detectors  for  the  electromagnet  LFI 

MRI  at  low  magnetic  field  (<10  mT)  without  cryogenic  or  hyperpolarization  techniques  presents 
unique  engineering  challenges.  Imaging  coils  must  maximize  coverage  over  the  volume  of  inter¬ 
est  while  minimizing  losses  in  a  regime  unusual  in  contemporary  MRI — where  Johnson  noise 
dominates  the  noise  floor.  Our  previous  work  validated  our  scanner  hardware  and  pulse  se¬ 
quences  MRI  at  6.5  mT.  The  purpose  of  the  present  work  is  to  construct  a  high-performance  coil 
for  low- field  imaging  of  the  human  brain  in  vivo. 


Spiral  volume  helmet  detectors:  receive  coil  development 


Our  previous  coils  focused  on  implementing  array  coil  technology  common  at  clinical  field 
strengths,  and  we  constructed  an  8  channel  receive-only  coil  for  our  6.5  mT  low-field  imager  at 
276  kHz.  While  our  results  were  promising  [11],  [12],  a  lack  of 
low  noise  high-impedance  pre-amplifiers  at  our  frequency  re¬ 
gime  prevented  us  from  achieving  sufficient  coil  decoupling, 
resulting  in  images  with  poor  SNR.  We  present  here  a  single 
channel  coil  with  greatly  improved  SNR  for  these  applications. 


We  designed  a  spiral  volume  coil  design  for  combined  Tx  and 
Rx  and  optimized  for  human  head  imaging  in  the  transverse  BO 
field  of  our  low  field  scanner.  This  places  the  receive  element 
close  to  the  head  while  still  covering  the  entire  volume  uniform¬ 
ly.  Additionally,  the  symmetry  of  a  spherical  spiral  results  in  a 
homogeneous  magnetic  field  over  its  volume,  minimizing  B1 
inhomogeneity  [13],  [14],  Losses  in  low  field  imaging  are  mainly 
due  to  resistive  losses  in  the  coil  (the  so-called  Johnson  noise  regime).  Rather  than  lower  this  re¬ 
sistance  by  cooling,  we  instead  use  of  multistranded  Litz  wire.  When  compared  to  a  solid  copper 
wire  of  the  same  gauge,  a  loop  made  of  Litz  wire  will  have  the  same  inductance  but  a  fraction  of 
the  AC  resistance,  resulting  in  a  lower  noise  floor  and  higher  image  SNR. 


Figure  8:  Form  fitting 
helmet  with  30-turn  spi¬ 
ral  design. 
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The  optimal  number  of  turns  and  appropriate  wire  gauge  were  determined  through  simulation. 
The  helmet  shape  was  designed  to  fit  closely  to  an  anatomically  accurate  human  head  model.  The 

spiral  path  was  generated  in  a 
custom  MATLAB  (Mathworks 
Inc.)  script  and  imported  into 
BiotSavart  (Ripplon  Software 
Inc.)  for  an  estimate  of  induct¬ 
ance.  Resistance  was  computed 
from  the  estimated  wire  length 
and  the  characteristic  imped¬ 
ance  per  length  of  wire  for  that 
gauge.  Coil  Q  was  estimated  as 
wL/R.  Bandwidth  was  estimated 
from  the  Q  value.  We  iterated 
through  several  designs  until  we 
achieved  our  target  Q~40  and 

BW=6  kHz. 

From  simulations,  we  built  a  30-turn  spiral  wound  with  Type  1  40/38  Litz  wire.  Our  aligned 
turn-to-turn  spacing  was  5.6  mm.  The  final  design  was  3D  printed  using  fused  deposition  model¬ 
ing  technology  with  polycarbonate  in  a  Fortus  360  me  printer  (Stratasys,  Eden  Prairie,  MN, 
USA),  (Figure  8). 

Images  were  acquired  at  6.5  mT  in  a  head  shaped  phantom  using  3D  b-SSFP  with  50%  under¬ 
sampling.  Resulting  voxel  size  was  3.9x3.5x16.7  mm3  (NA=110)  (Figure  9a).  A  maximum  in¬ 
plane  SNR  of  130  was  calculated.  An  image  of  a  similar  slice  obtained  with  the  8-channel  array 
coil  is  shown  for  comparison  (Figure  9b,  combined  image  using  RSS,  voxel  size  =  4.4x4.2x22 
mm3,  NA  =  80,  max  SNR  =  16.5). 

18  slices  of  a  head- shaped  structured  resolution  phantom  imaged  with  optimized  spiral  volume 
coil  are  shown  in  Figure  10.  This  four-quadrant  phantom  is  filled  will  liquid-filled  spheres  of  var¬ 
ious  diameters.  Voxel  size  in  this  scan  is  3. 5x2. 5x5  mm3  with  NA=40  and  maximum  SNR  of  24.8. 
Total  acquisition  time  is  10  minutes. 


Figure  9:  Equivalent  slices  in  a  head  phantom  acquired 
using  the  single  channel  spiral  coil  (A)  and  the  8-channel 
array  (B)  at  6.5  mT.  Images  are  scaled  by  their  maximum 
intensity. 
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Figure  10:  Axial  view  of  a  multi-compartment  head-shaped  phantom  acquired  with  the  Litz 
spiral  volume  coil  at  6.5  mT  in  10  minutes. 


This  30-turn  single  channel  volume  spiral  coil  designed  for  low  field  imaging  significantly  out¬ 
performed  our  8 -channel  array  coil.  By  minimizing  losses  in  the  coil,  maximizing  filling  factor 
and  eliminating  coupling  issues,  we  were  able  to  obtain  significantly  higher  SNR.  We  have  suc¬ 
cessfully  designed  and  constructed  a  high  filling  factor  human-head  Tx/Rx  spiral  volume  imag¬ 
ing  coil  for  6.5  mT  (276  kHz). 

OMRI  Injury  Imaging 

TASK  2A:  OMRI  Hardware  Development 

Overhauser- enhanced  MRI  (OMRI)  is  an  electron-proton  double  resonance  imaging  technique 
of  much  interest  due  to  its  ability  to  detect  the  concentration  and  distribution  of  free  radicals. 
Tracking  of  exogenous  free  radicals  with  OMRI  in  vivo  has  enabled  the  development  of  oxymetry 
probes  [15]  and  the  imaging  of  redox  reactions  [16].  The  large  gyromagnetic  ratio  of  electrons 
(28  GHz/T)  demands  that  in  vivo  OMRI  is  performed  at  very  low  magnetic  fields  (~10  mT)  in 
order  to  minimize  RF  heating  and  penetration  depth  issues.  Operation  at  low  magnetic  field 
causes  a  drastic  reduction  in  NMR  sensitivity  despite  the  signal  enhancement  that  comes  from 
the  Overhauser  effect,  and  emphasizes  the  need  for  high  S/N  probes.  OMRI  probe  design  is  still 
relatively  unexplored,  despite  its  importance,  and  presents  challenges  unique  to  the  frequencies 
of  operation  (fH  =  276  kHz  and  fe  =  140.8  MHz  in  our  experiments  at  6.5  mT). 

A  key  challenge  in  implementing  OMRI  for  TBI  will  be  detecting  the  small  concentrations  of 
radical  present  in  vivo.  We  have  previously  shown  that  10  p(M  nitroxide  radical  in  water  and  50 
piM  nitroxide  radical  in  water  can  be  detected  with  Overhauser  spectroscopy  (Figure  11)  and 
OMRI  (Figure  12),  respectively.  The  concentration  of  free  radical  at  the  location  of  injury  re- 
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mains  unknown ,  but  estimates  range  from  50  nM  to  50  /uM  [17] -[24].  Additionally,  detecting  free 
radicals  in  vivo  will  be  more  difficult  than  in  phantoms  due  to  extra  relaxation  pathways  reducing 
the  efficiency  of  DNP.  Therefore  we  would  like  to  try  and  improve  our  detection  efficiency  even 
further. 


Figure  11:  Measurement  of  the  sensitivity  threshold  to  small  concentrations  of  nitroxide  radical 
via  Overhauser  enhanced  NMR  spectroscopy.  These  spectra  were  acquired  in  the  electromagnet 
LFI  at  6.5  mT  (276  kHz).  The  ESR  frequency  was  140  MHz. 


Figure  12:  Photograph  of  phantom  used  to  demonstrate  the  sensitivity  of  OMRI  to  nitroxide  radi¬ 
cal  concentration.  All  seven  vials  have  very  similar  image  magnitudes  in  the  conventional  MRI 
image.  The  OMRI  scan  demonstrates  marked  image-based  free  radical  sensitivity.  The  signal 
magnitude  from  the  control  vial  remains  unchanged  in  OMRI  scan.  OMRI  enhancement  image 
is  computed  from  the  ratio  of  OMRI  to  MRI  magnitude.  MRI  and  OMRI  were  both  obtained  at 
6.5  mT.  The  NMR  frequency  is  276  kHz  and  the  ESR  frequency  is  140  MHz. 


Improving  radical  detection  sensitivity  through  double  resonator  design 

We  report  here  the  development  of  a  high  performance  OMRI  probe  with  broad  tunability,  ca¬ 
pable  of  imaging  enhancement  over  a  wide  range  of  free  radical  electron  ^-factors.  This  OMRI 
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probe  was  developed  because  it  was  determined  that  significant  improvement  in  our  OMRI 
hardware  were  needed  before  we  could  attain  sensitivity  to  superoxide.  We  have  developed  an 
OMRI  probe  consisting  of  an  NMR  solenoid  inside  a  modified  Alderman-Grant  Resonator 
(Figure  13),  to  interfaces  with  the  6.5  mT  electromagnet  LFI  (fH  =  276  kHz).  Images  were  ac¬ 
quired  using  our  recently  developed,  fast,  high-resolution  b-SSFP  based  OMRI  methodology 
[25].  NMR  probe  design  in  the  low-field  regime  is  fundamentally  different  to  that  at  conventional 
MRI  fields  as  thermal  noise  due  to  the  intrinsic  resistance  of  the  pickup  coil  dominates  over  sam¬ 
ple  noise.  This  leads  to  a  compromise  where  S/N  improvements  come  at  the  expense  of  imaging 
bandwidth  S/N  ~1/VBW.  An  85  turn  solenoid  was  wound,  using  low  AC  resistance 
5/39/42  litz  wire,  on  a  3D  printed  polycarbonate  former.  This  high  filling  factor  coil  has  a  band¬ 
width  of  3  kHz. 
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Figure  13:  NMR/ESR  rat  head  probe  for  operation  at  6.5  mT.  The  ESR  resonator  is  tuned  to  140.8 


MHz.  The  litz  wire  NMR  solenoid  coil  (276  kHz)  resonator  board  is  not  shown. 


A  modified  Alderman-Grant  ESR  resonator  was  built  using  copper  foil  on  Pyrex  tubing.  All  met¬ 
al  placed  in  close  proximity  to  the  NMR  solenoid  strongly  couples,  reducing  the  NMR  sensitivity. 
We  therefore  minimized  amount  of  copper  in  the  ESR  resonator.  Windows  were  removed  from 
the  panels  on  the  sides  of  an  Alderman- Grant  resonator,  a  region  of  low  current  flow  [26],  in  an 
attempt  to  reduce  coupling  whilst  maintaining  Bi  homogeneity.  Shielding  at  the  ends  of  the  reso¬ 
nator  prevents  high  electric  fields  at  the  capacitors  penetrating  the  imaging  volume,  important 
because  PRF-absorbed  ~  E2.  Slits  in  the  shielding  prevent  the  formation  of  closed  loops  that  couple  to 
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the  solenoid.  Figure  14  demonstrates  the  high  Bi  homogeneity  and  strong  E  suppression  in  our 
ESR  resonator. 


Figure  14:  COMSOL  Multiphysics  simulations  (a)  Our  modified  Alderman-Grant  Resonator,  (b) 
Plot  of  |E|  showing  that  the  electric  field  is  strongly  suppressed  inside  the  resonator,  (c)  Bie  shows 
high  homogeneity  in  the  imaging  volume  -  less  than  10%  variation  across  the  imaging  region. 


TEMPOL  (4-hydroxy-TEMPO)  a  stable  radical,  is  detected  by  OMRI  with  high  sensitivity,  and 
because  of  its  small  size  may  facilitate  imaging  blood  brain  barrier  leakage  in  cases  of  oxidative 
stress  [27].  Simulations  were  validated  in  2  mM  TEMPOL  solutions  using  simple  spectroscopic 
measurements  as  well  as  fast  imaging  strategies  [25].  The  ESR  resonator  was  tuned  to  141  MHz, 
the  lowest  frequency  of  the  TEMPOL  triplet  state  at  6.5  mT,  to  minimize  RF  heating  during  ESR 
irradiation. 


In  vivo  performance 

Previously  we  described  the  development  of  an  OMRI  probe  for  a  rat  head  model  utilizing  an 
NMR  solenoid  and  ESR  surface  coil.  Testing  of  this  new  OMRI  probe  indicates  that  it  has  3x  the 
S/N  of  the  probe  reported  [27],  and  rectifies  problems  with  Bie  homogeneity,  yielding  homoge¬ 
neous  enhancement  of  -6.7  in  2  mM  TEMPO  when  10  W  of  RF  power  is  applied.  In  vivo  OMRI 
signal  enhancement  is  clearly  visible  in  the  rat  brain  after  TEMPOL  injection,  as  shown  in  Figure 
15.  As  the  Overhauser-enhanced  signal  has  a  phase  opposite  to  that  of  the  thermal  signal,  the 
phase  image  in  Figure  15  provides  sensitive  contrast  in  regions  of  low  TEMPOL  concentration. 
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Figure  15:  OMRI  images  acquired  from  a  rat  at  6.5  mT  following  injection  of  1  mL  of  150  mM 
TEMPOL.  Five  slices  from  an  11  slice  data  set  shown.  OMRI  (NAM)  imaging  time  was  9  se¬ 
conds.  Anatomical  MRI  (NA=30)  was  acquired  in  the  OMRI  scanner  with  ESR  power  disabled. 

All  images,  voxel  size:  1.1  x  1.6  x  8  mm3,  Matrix:  128  x  35  x  11.  A  fiducial  containing  2  mM  of 
TEMPOL  was  placed  outside  the  animal  and  is  seen  in  both  MRI  and  OMRI  images. 

In  summary,  we  have  implemented  a  high  performance  probe  for  high  temporal  and  spatial  reso¬ 
lution  OMRI  and  demonstrated  it  in  vivo  in  the  brain  of  living  rats.  The  S/N  of  this  probe  maybe 
further  improved,  whilst  maintaining  imaging  bandwidth,  by  using  an  active  feedback  circuit 
[28].  Higher  S/N  could  also  be  realized  through  supercooling  of  the  NMR  solenoid  or  by  using  a 
free  radical  with  a  longer  Tie  such  as  triphenylmethyl  [15].  This  probe  may  allow  the  in  vivo  de¬ 
tection  of  rapid  redox  changes  in  pathologic  tissues,  specifically  in  the  context  of  brain  trauma  or 
stroke. 


KEY  RESEARCH  ACCOMPLISHMENTS: 

Electromagnet  6.5  mT LFI: 

•  Designed,  constructed,  and  characterized  optimized  high-efficiency  single  channel  spiral 
“helmet”  for  human  head  imaging  at  6.5  mT 

•  Designed  and  built  new  multi- compartment  head-shaped  imaging  phantom 

•  Acquired  high-quality  images  from  multi-compartment  phantom  using  human  head  coil 
2. 5x3. 5x5  mm  resolution  (18  slices)  in  10  minutes 

Permanent  magnet  Halbach  LFI: 

•  Simulated  TRASE  array  to  enable  3D  encoding 
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•  Designed  and  built  TRASE  array  coils  with  switchable  180°  phase  shifter 

•  Implemented  3D  encoding  using  simultaneous  2D  rSEM  and  TRASE 

•  Implemented  WURST-40  acquisition  to  mitigate  quadratic  phase  from  frequency- swept 
pulses 

•  Demonstrated  proof-of-concept  3D  imaging  in  the  inhomogenous  45  kilogram  Halbach 
magnet  without  the  use  of  gradient  coils. 

Overhauser  MRI  in  electromagnet  LFI: 

•  Established  collaboration  with  D.  Whitney  King  (Colby  College)  to  provide  superoxide 
generator  for  OMRI  sensitivity  calibration 

•  Developed  and  built  high  performance  OMRI/MRI  probe  with  wide-range  OMRI  tuning 
and  low  SAR,  and  3-fold  improved  MRI  SNR. 

•  Acquired  11  slice,  l.lxl. 6x8  mm  in-plane  resolution  MRI  (4.5  minutes)  and  OMRI  (9 
seconds)  in  rats 

•  Developed  methodology  to  measure  dynamic  time  resolved  OMRI  in  vivo 

REPORTABLE  OUTCOMES; 

1.  C  LaPierre,  M  Sarracanie,  B  D  Armstrong,  J  Polimeni,  and  M  S  Rosen,  Overhauser- 
enhanced  MRI  with  SENSE  Acceleration  in  the  Johnson  Noise  Dominated  Regime,  pre¬ 
sented  at  the  55th  Meeting  of  the  Experimental  NMR  Conference,  Boston,  MA,  March  23- 
28  2014. 

2.  M  Sarracanie,  B  D  Armstrong,  and  M  S  Rosen,  High  Speed  MR  Fingerprinting  at  6.5  mT, 
presented  at  the  55th  Meeting  of  the  Experimental  NMR  Conference,  Boston,  MA,  March 
23-28  2014. 

3.  M  Sarracanie,  P  Stanwix,  and  M  S  Rosen,  A  Single  Channel  Variable  Density  Spiral  Coil 
for  High  Sensitivity  MRI  at  6.5  mT,  presented  at  the  55th  Meeting  of  the  Experimental 
NMR  Conference,  Boston,  MA,  March  23-28  2014. 

4.  C  Z  Cooley,  J  P  Stockmann,  M  Sarracanie,  B  D  Armstrong,  M  S  Rosen,  and  L  L  Wald,  2D 
Imaging  in  a  Portable  MRI  Scanner  without  Gradient  Coils,  presented  at  the  55th  Meeting 
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CONCLUSION; 


The  reason  MRI  is  not  widely  deployable  is  that  high-strength  magnetic  fields  (of  order  1  T)  are 
necessary  with  conventional  MRI  to  obtain  useful  brain  images.  Such  high  field  scanners  involve 
large,  heavy,  fragile,  expensive  equipment  (such  as  superconducting  magnets)  that  are  difficult  to 
site  in  field  hospitals.  We  contend  that  that  low-magnetic-field  implementations  of  MRI  can  be 
developed  to  allow  robust,  transportable  imaging  modalities  well  suited  to  diagnose  the  types  of 
battlefield  injuries  prevalent  in  TBI  and  practical  for  operation  in  field  hospitals.  Application  of 
the  suite  of  techniques  and  technologies  from  our  work  could  advise  future  development  of  a  de¬ 
ployable  device  with  a  high  diagnostic  impact  and  could  be  transformative,  enabling  improved 
diagnosis  and  monitoring  of  battlefield  injuries  prevalent  in  TBI. 

The  electromagnet  LFI  provides  an  ideal  state-of-the  art  test  bed  for  all  of  the  novel  acquisition, 
detection  methodologies,  and  reconstruction  algorithms  including  navigators  and  sparse  sam¬ 
pling,  and  additionally  will  provide  necessary  experience  and  data  to  advise  optimal  construction 
and  magnetic  field  for  any  future  electromagnet-based  deployable  systems.  In  its  current  configu¬ 
ration,  this  state-of-the-art  scanner  enables  high-performance  spectroscopy  and  8-channel  imag¬ 
ing  at  6.5  mT,  and  is  fully  equipped  for  Overhauser  DNP  experiments.  We  continue  to  improve 
the  ability  to  acquire  images  based  on  the  intrinsic  'H  signal  by  combining  new  low-field  hard¬ 
ware  methodology  and  advanced  pulse  sequence  and  reconstruction  methods. 

Our  second  low  field  imaging  test  bed,  the  permanent  magnet  Halbach  LFI  is  based  around  a 
lightweight  (45  kg)  and  portable  magnet  with  a  built-in  encoding  field  for  MRI.  It  is  a  highly  spe¬ 
cialized  and  potentially  disruptive  scanner  that  could  greatly  ease  both  the  cost  and  burden  of  a 
field-forward  instrument  purpose-built  for  TBI  imaging.  The  2D  imaging  results  acquired  in  the 
highly  inhomogeneous  magnetic  field  of  the  Halbach  scanner  using  simultaneous  NMR  meas¬ 
urement  to  track  magnet  drift,  and  multi-channel  receive  arrays  to  unwrap  image  aliasing,  and 
are  encouraging.  A  truly  portable  MR  system  based  on  this  technology  has  the  potential  to  quick¬ 
ly  detect  brain  injury  at  the  site  of  injury,  for  instance  in  an  ambulance  prior  to  transportation  to 
the  hospital. 

We  have  also  been  developing  technology  for  a  wholly  new  approach  to  brain  imaging  that  may 
revolutionize  the  use  of  MRI  for  the  assessment  and  treatment  of  secondary  brain  injury  follow¬ 
ing  TBI:  direct  tomographic  detection  of  endogenous  free  radicals  as  an  early  marker  for  TBI. 
This  new  form  of  contrast,  “injury- sensitive  MRI”,  is  based  on  the  enormous  signal  enhancement 
attainable  with  Overhauser  DNP,  converting  the  electron  spin  of  endogenous  free  radicals  into 
nuclear  polarization  using  the  Overhauser  effect  and  subsequently  imaging  that  modified  nuclear 
polarization  using  low- field  MRI  (OMRI). 
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The  imaging  experiments  presented  above  in  the  optimized  electromagnet  LFI  are  compelling: 
high  quality  'H  brain  MRI  at  6.5  mT  is  in  fact  attainable  over  reasonable  averaging  times  using 
high  performance  hardware  and  advanced  MRI  sequences.  A  critical  question  for  this  work  is  the 
understanding  of  the  clinical  balance  and  impact  that  the  tradeoff  between  acquisition  time  and 
resolution  implicit  in  imaging  plays  on  the  operation  of  a  deployable  scanner. 

The  low-field  OMRI  results  presented  above  demonstrate  promise  for  high-speed  high- 
resolution  free-radical  imaging,  and  offers  new  perspectives  for  the  measurement  of  free  radicals 
in  living  organisms.  The  long-term  goal  of  this  work  is  an  in  vivo  implementation  of  this  tech¬ 
nique.  Free-radical  sensitive  low-field  OMRI  as  a  novel  imaging  and  diagnostic  MRI-based 
method  with  specificity  to  secondary- injury  has  the  potential  to  clarify  the  mechanisms  involved 
in  secondary  damage  and  the  local  effects  of  novel  therapies.  Time-critical  imaging  of  free  radi¬ 
cals  can  provide  the  unique  interventional  access  critically  needed  for  drug  therapies  and  this 
technique  would  fill  the  clear  need  for  an  “injury  imager”,  suitable  for  non-invasive  tomographic 
measurement  (Figure  16). 


Figure  16:  Schematic  of  the  development  of  a  new  interventional  imaging  and  treatment  tool  for 
in  vivo  detection  of  endogenous  free  radicals  associated  with  secondary  injury  following  TBI  us¬ 
ing  low-field  OMRI  as  developed  here.  This  same  approach  can  similarly  open  up  new  treatment 
windows  into  stroke  and  other  ischemic  events.  MRI-based  in  vivo  free  radical  imaging  using 
OMRI  is  impossible  at  high-field  due  to  the  inability  of  the  ESR  pulse  to  penetrate  into  tissue,  and 
the  danger  of  RF  heating  in  the  microwave  regime. 

Successful  demonstration  of  this  system  in  vivo  will  provide  image-based  injury  specificity  for 
TBI  and  for  the  first  time  open  a  window  into  a  critical  mechanism  of  TBI-related  disease  and  an 
unambiguous  non-invasive  in  vivo  marker  for  cerebral  injury.  Time-critical  imaging  of  free  radi¬ 
cals  associated  with  secondary  imaging  can  provide  the  unique  interventional  access  critically 
needed  to  develop  and  deliver  drug  therapies  to  mitigate  brain  damage,  reducing  disability  and 
death  from  secondary  injury  following  TBI.  The  low-field  approach  would  allow  placement  of 
this  critical  new  diagnostic  “injury  imager”  tool  close  to  the  source  of  injury — as  well  as  in  hospi¬ 
tals  and  clinics  where  rehabilitation  takes  place. 
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APPENDICES: 


Year  4  Statement  of  Work: 

We  attach  for  reference  the  SOW  for  the  No  Cost  Extension  of  Y4  as  Appendix  1. 

Abstracts  accepted  for  presentation: 

We  attach  seven  abstracts  (Appendices  2-8)  accepted  for  presentation  in  2015  that  result  from 
research  presently  funded  in  the  Rosen  lab  by  DoD/DMRDP.  Abstracts  marked  with  a  $  will  be 
presented  as  a  poster  at  the  56th  Experimental  NMR  Conference  (ENC)  to  take  place  20-24  May 
2015  in  Pacific  Grove,  CA,  and  abstracts  marked  with  a  §  will  be  presented  as  a  poster  at  the  24th 
Annual  Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine  (ISMRM)  to 
take  place  1-5  June  2015  in  Toronto,  Canada.  A  symbol  printed  in  red  indicates  that  the  abstract 
will  be  presented  as  a  talk. 

1.  3D  Balanced-EPI  Magnetic  Resonance  Fingerprinting  at  6.5  mT 

2.  Low  Cost  High  Performance  MRI 

3.  $  §  Dynamic  in  vivo  free  radical  imaging  with  Overhauser-enhanced  MRI 

4.  3D  Imaging  in  a  Portable  MRI  Scanner  using  Rotating  Spatial  Encoding  Magnetic 
Fields  and  Transmit  Array  Spatial  Encoding  (TRASE) 

5.  $  §  High  Performance  Probe  for  in  vivo  Overhauser  MRI 

6.  A  single  channel  spiral  volume  coil  for  in  vivo  imaging  of  the  whole  human  brain  at 
6.5  mT 

7.  §  Transmit  Array  Spatial  Encoding  (TRASE)  with  broadband  WURST  pulses  for  robust 
spatial  encoding  in  inhomogeneous  BO  fields 

Manuscripts  published/accepted  for  publication: 

We  attach  as  Appendix  9  our  manuscript  “2D  imaging  in  a  lightweight  portable  MRI  scanner 
without  gradient  coils”,  published  in  Magnetic  Resonance  in  Medicine  in  2014. 
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HARVARD  UNIVERSITY  DEPARTMENT  OF  PHYSICS 

1 7  Oxford  Street 

CAMBRIDGE,  MASSACHUSETTS  02138 


13  December  2013 

Lance  L.  Nowell 
Contract  Specialist 
Assistance  Branch  4 

U.S.  Army  Medical  Research  Acquisition  Activity 

Attn:  MCMR-AAA-AD 

820  Chandler  Street 

Fort  Detrick,  MD  21702-5014 

RE:  NCE  Request  for  ARO  DMRDP  Award  #  W81XWH-11-2-0076 


Per  you  email  request  of  Tues  December  12  please  find  the  following  information  in  support  of  my 
EWOF  request. 

During  the  extension  period,  we  will  be  completing  the  following  tasks: 

*  Part  of  Task  1A:.  Enable  3rd-axis  (axial)  encoding  with  Halbach  permanent  imager  using  TRASE 

coil.  This  is  a  critical  piece  of  technology  needed  to  fully  realize  true  3D  MRI  in  a  portable  lightweight 
(<100  lb)  head  imager,  and  is  one  of  the  two  novel  MRI  scanners  developed  under  this  grant.  This  task 
been  delayed  because  of  technical  issues  surrounding  RF  power  handling  issues. 

*  Part  of  task  1C:  3.  Test  image-based  free  radical  sensitivity  of  OMRI  system.  We  have  identified 
Professor  of  Chemstry  D.  Whitney  King  (Colby  College)  as  a  collaborator  on  this  task.  His  laboratory  has 
developed  a  calibrated  superoxide  generator,  which  he  will  bring  to  the  Rosen  Laboratory  at  Harvard.  We 
will  use  this  instrument  to  directly  investigate  OMRI  sensitivity  to  biologically  relevant  superoxide  at  na¬ 
nomolar  to  micromolar  concentrations.  This  is  a  critical  step  to  understand  the  biological  relevance  and 
applicability  of  the  new  method  of  free  radical  imaging  developed  in  the  Rosen  lab  under  this  grant.  This 
task  had  been  delayed  because  more  time  then  budgeted  was  spent  improving  detection  strategies,  as  we 
continued  to  make  gains  in  this  regime. 

Plan  for  use  of  unobligated  funds: 

This  no-cost  extension  is  requested  because  of  technical  delays  in  the  current  program.  This  request  for  a 
no  cost  extension  will  enable  the  previously  approved  objectives  to  be  carried  out.  During  the  requested 
extension,  the  budget  will  be  spent  in  accordance  with  the  original  budget  and  budget  justification.  No 
additional  funds  are  requested. 

With  very  best  regards, 


Matthew  S.  Rosen,  Ph.D. 


LOW-FIELD  MRI  AND  H YPERPOLARIZED  MEDIA  LABORATORY 
MGH/M  ARTI  NOS  CENTER  FOR  BIOMEDICAL  IMAGING 
149  1  3TH  STREET,  SUITE  2301 
CHARLESTOWN,  MA  02  1  29 
LAB:  617  643  8636,  FAX:  617-726-7993 
martinos . org/lf i 
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3D  Balanced-EPI  Magnetic  Resonance  Fingerprinting  at  6.5  mT 

Mathieu  Sarracanie1’2,  Ouri  Cohen1,  and  Matthew  S  Rosen1,2 

1MGH/A.A.  Martinos  Center  for  Biomedical  Imaging,  Charlestown,  MA,  United  States,  2Department  of  Physics,  Harvard  University,  Cambridge,  MA,  United  States 

Purpose 

In  recent  work  [1],  we  demonstrated  high  speed  MRI  in  the  very  low  magnetic  field  regime  (6.5  mT)  using  a  balanced  steady  state  based  (b-SSFP 
[2])  acquisition  scheme.  b-SSFP  provides  the  highest  SNR  per  unit  time  [2]  and  image  contrast  depends  on  the  ratio  T/iy  At  very  low  field,  most 
species  have  T2  relaxation  times  approaching  Tb  so  b-SSFP  images  are  essentially  proton  density  (PD)  weighted.  In  previous  work  [3]  we  have 
shown  that  2D  MR  Fingerprinting  [4]  can  be  implemented  at  low  magnetic  field  and  provide  simultaneous  quantification  of  Ti  and  T2  as  well  as  pro¬ 
ton  density  and  B0  field  maps.  MRF  at  low  magnetic  field  creates  a  rapid  dynamic  series  of  low  signal  to  noise  ratio  (SNR)  images  where  the  magni¬ 
tude  of  each  voxel  of  each  image  changes  at  every  time  step.  Generally,  the  TR  and  flip  angle  of  each  image  in  the  time  series  is  varied  pseudo- 
randomly  [5].  Here,  we  demonstrate  MRF  in  3D  at  6.5  mT,  using  an  optimized  set  of  15  flip  angles  and  repetition  times  (FA/TR),  in  a  Cartesian  ac¬ 
quisition  of  &-space  with  a  new  hybrid  b-SSFP-EPI  sequence. 

Methods 

The  low  field  MRI  scanner  was  previously  described  [5].  The  imaging  sequence  is  a  hybrid  multishot  b-SSFP-EPI  with  an  echo-train  of  4  echoes 
tr  (Fig.  1).  The  sequence  was  set  with  matrix  size  =  64x64x5,  corresponding  voxel  size  = 

t  (2.5x3.5x10)  mm3,  FOV  =  (125x175x50)  mm3,  number  of  average  NA  =  2.  The  minimum  TR  was 

62  ms  with  9091  Hz  bandwidth.  The  total  acquisition  time  was  25  min.  A  flip  (FA/TR)  trajectory  of 
length  N=15  was  generated  using  an  optimization  method  previously  described  [6].  Lack  of  SNR  at 
low  magnetic  field  required  redesigning  our  optimization  scheme  with  a  narrow  range  of  larger  flip 
angles  [60-120°]  and  shorter  TRs  [62-400  ms].  The  optmized  FA/TR  trajectory  is  compared  to  a 
non-optimized  trajec- 


Figure  1:  Pulse  sequence  diagram  of  the  hybrid 
b-SSFP-EPI  sequence  used  for  3D  MR  Finger¬ 
printing.  The  FAs  and  TRs  were  set  according  to 
the  optimized  trajectory. 


tory  in  Fig.  2.  The 
imaged  phantom  con¬ 
sists  in  a  stack  of 
three  10  mm  thick 
compartments  of 
similar  volume  with 
different  T,  and  T2 
properties. 
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Results 

Each  image  generated  in  the 
reconstructed  fingerprinting 
set  (Figure  3)  reveals  differ¬ 
ent  information.  In  regions 
with  no  signal  (Slice  1,  Slice 
5),  no  matched  value  can  be 
found.  The  spin  density  (M0)  map  of  Figure  3. a  is  equivalent  to  traditional  b-SSFP,  and  only  subtle  differ¬ 
ences  are  seen  between  compartments.  Figure  3.  b-c  reveals  that  the  compartment  in  slice  3  has  the  lowest 
T |  and  T2  values,  and  that  the  compartment  in  slice  4  has  the  highest  T,  and  T2  values.  Mean  T!  values  are 
200ms,  168ms,  and  320  ms  in  slices  2,  3,  and  4  respectively.  Mean  T2  values  are  180,  157,  and  260  ms  in 
slices  2,  3,  and  4  respectively. 


Figure  2:  Comparison  of  our  optimized 
trajectory  with  a  random  non-optimized 
trajectories  for  N=15. 


Figure  3  :  MR  Fingerprinting  results  at  6.5  mT  :  a.  M0,  b.  Ti,  and  c.  T2  in  3D,  in 
a  3  compartment  phantom.  Each  slice  in  the  figure  matches  one  of  the  phantom 
compartments. 


Conclusion 

We  have  demonstrated  3D 
MR  Fingerprinting  at  very 
low  magnetic  field  with  a 

hybrid  b-SSFP-EPI  sequence  enabling  fast  and  robust  acquisition  of  k-space.  The  optimized  FA/TR  strate¬ 
gy  provides  good  dispersion  while  drastically  reducing  the  total  acquisition  time.  We  measure  quantitative 
parameters  in  3D,  and  generate  several  image  contrasts  in  a  single  acquisition  (proton  density,  Tls  T2)  in  less  than  30  minutes.  This  technique  is  of 
particular  relevance  at  low  magnetic  field  where  SNR  and  contrast  are  tied  to  long  acquisition  times.  The  combination  of  3D  MRF  with  low  field 
MRI  scanners  has  great  potential  to  provide  clinically  relevant  contrast  with  portable  low  cost  MR  scanners. 

References:  [1]  Sarracanie  M  et  al.  Proc.  ISMRM  2013  #5322;  [2]  Scheffler  K  et  al.  Eur  Radiol  2003  13:2409-18;  [3]  Sarracanie  M  et  al.  Proc. 
ISMRM  2014  #6370;  [4]  Ma  D  et  al  Nature  2013  495:187-193;  [5]  Tsai  LL  et  al.  JMR  2008;  193:174-85;  [6]  Cohen  et  al,  Proc  ISMRM  2014  #7153 

Ackowledgement:  This  research  was  supported  by  the  Department  of  Defense,  Defense  Medical  Research  and  Development  Program,  Applied  Re¬ 
search  and  Advanced  Technology  Development  Award  W81XWH-1 1-2-0076  (DM09094). 
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Low  Cost  High  Performance  MRI 

Mathieu  Sarracanie1,2,  Cristen  LaPierre1,2,  Najat  Salameh1,2,  David  E  J  Waddington1’3,  Thomas  Witzel1,  and  Matthew  S  Rosen1,2 
'  MGH/A.A.  Martinos  Center  for  Biomedical  Imaging,  Charlestown,  MA,  United  States,  2Department  of  Physics,  Harvard  University,  Cambridge,  MA,  United  States, 
3 ARC  Centre  of  Excellence  for  Engineered  Quantum  Systems,  School  of  Physics,  University  of  Sydney,  Sydney,  NSW,  Australia 


Purpose 

Magnetic  Resonance  Imaging  (MRI)  is  unparalleled  in  its  ability  to  non-invasively  visualize  anatomical  structure  and  function  with  high  spatial  and 
temporal  resolution  and  a  broad  range  of  clinically  relevant  contrasts.  Yet  to  overcome  the  low  sensitivity  inherent  in  NMR-based  detection,  the  vast 
majority  of  clinical  scanners  incorporate  superconducting  magnets  operating  at  1.5,  3  tesla  (T),  and  more  exceptionally  at  7T.  These  powerful 
magnets  are  massive,  costly  to  purchase  and  maintain,  and  operate  with  very  strict  infrastructure  demands  and  siting  requirements  that  preclude 
operation  in  many  environments.  For  brain  imaging,  low-field  alternatives  have  been  developed  that  rely  on  superconducting  quantum  interference 
devices  (SQUID)  detector  arrays  combined  with  magnetic  field  cycling  techniques  [1-4].  To  date,  this  approach  is 
limited  by  prohibitively  long  acquisition  times  (>  1  hr)  and  restrictions  on  attainable  fields  of  view.  We  demonstrate 
here  3D  brain  MRI  in  vivo  at  6.5mT  (more  than  450  times  lower  than  clinical  MRI  scanners)  in  6  minutes,  by 
combining  a  high  performance  single  channel  Tx/Rx  coil  with  modern  undersampling  techniques  and  b-SSFP  [5]. 

Methods 

To  minimize  B,  inhomogeneity  and  maximize  filling  factor,  a  close  fitting  single  channel  inductive  head  coil  for 
operation  at  276  kHz  was  designed  and  3D  printed  (Fig.  1).  The  coil  features  a  30-turn  3D  spiral  with  turn-to-turn 
distance  of  5.6  mm,  thus  ensuring  the  B]  field  produced  is  everywhere  orthogonal  to  the  main  magnetic  field  B0.  The 
hemispheric  spiral  design  results  in  a  very  homogeneous  magnetic  field  [6,7]  over  the  volume  of  interest,  making  it 
suitable  for  both  RF  transmit  and  receive.  Litz  wire  was  selected  due  to  its  lower  AC  resistance  compared  to  solid 
copper  wire,  reducing  the  Johnson  noise  of  the  coil  while  maintaining  the  same  inductance.  Previously  [8],  we 
described  our  use  of  undersampling  strategies  to  accelerate  low-field  imaging.  We  make  use  of  this  here  by  randomly 
sampling  50%  of  k-space  using  a  variable  density  Gaussian  pattern.  Once  reconstructed,  the  images  were  apodized 
and  processed  using  Perona  and  Malik  anisotropic  diffusion  filtering  [9]  (ADF). 

Results 

Three-dimensional  undersampled  images  acquired  in  6  minutes  (NA=30)  are  shown  in  Fig.  2  for  each  of  the  three  spatial  orientations  (axial,  coronal, 
and  sagittal).  We  obtain  2.5x3.5x8.5  mm3,  2.5x3.5x14.4  mm3,  and  2.5x3.5x11.5  mm3  voxel  resolution  images  in  axial,  sagittal  and  coronal 
orientation,  respectively.  Images  with  NA=160  acquired  at  6.5  mT  are  compared  to  b-SSFP  images  acquired  at  3T  in  Fig.  3. 


Figure  1  :  3D  printed  single 
channel  volume  head  coil. 
The  coil  was  used  for  both 
transmit  and  receive  at  276 
kHz  (6.5  mT) 


Conclusion 

At  6.5  mT,  (more  than 
450  times  lower  than 
clinical  MRI  scanners)  we 
demonstrate  2.5x3.5x8.5 
mm3  resolution  in  the 
living  human  brain  in  6 
minutes.  We  contend  that 
robust  non-field-cycled 
low-field  implementations 
of  MRI  (<  10  mT)  have 
the  potential  to  make 
clinically  relevant  images 
and  set  new  standards  for 
a  completely  new 
category  of  affordable  and 
robust  portable  devices. 
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Figure  2:  3D  images  of  the  living  brain  (7  central  slices  of  15  slice  dataset  are  shown) 
acquired  in  6  minutes  at  6.5  mT  (276  kHz)  in  a.  axial,  b.  coronal,  and  c.  sagittal 
orientation.  Acquisition  matrix:  64x75x15,  Flip  angle  =  70°,  BW  =  9091  Hz,  TR/TE= 
23/11.7  ms,  NA=30,  voxel  size:  a.  (2. 5x3. 5x8. 5)  mm3,  b.  (2.5x3.5x11.5)  mm3,  and  c. 
(2.5x3.5x14.4)  mm3. 


Figure  3:  Comparison  of  b-SFFP  images 
at  a.  3T  and  b.  6.5  mT.  Strong  banding 
artifacts  appear  at  high  magnetic  field 
(yellow  arrows)  whereas  no  artifact  is  seen 
in  the  images  acquired  at  low  field. 


References:  [1]  Zotev  et 
al.  J.  Magn.  Res.  2008 

194(1):  1 15-120;  [2]  Espy  et  al  J.  Magn.  Res.  2012  228:1-15;  [3]  Yesanen  et  al  Magn.  Reson.  Med.  2012  69(6):1795-1804;  [4]  Inglis  et  al  PNAS 
2013  110(48):  19194-201;  [5]  Scheffler  K  et  al  Eur  Radiol  2003  13:2409-18;  [6]  Harpen  J.  Magn.  Res.  1991  94(3):550-556;  [7]  Everett  et  al  J. 
Scientific  Instruments  1966  43:470-474;  [8]  Sarracanie  et  al  Magn.  Reson.  Med.  2013  71:735-745;  [9]  Perona  et  al  1990  IEEE  Trans.  Pattern  Anal. 
Mach.  Intel!  12(7):629-639. 
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Dynamic  in  vivo  free  radical  imaging  with  Overhauser- enhanced  MRI 

Mathieu  Sarracanie1’2,  Fanny  Herisson3,  Najat  Salameh1,2,  David  E  J  Waddington1,4,  Cenk  Ayata3,  and  Matthew  S  Rosen1,2 
'MGH/A.A.  Martinos  Center  for  Biomedical  Imaging,  Charlestown,  MA,  United  States,  2Department  of  Physics,  Harvard  University,  Cambridge,  MA,  United  States, 

3 Neurovascular  Research  Lab,  Department  of  Radiology,  Massachusetts  General  Hospital,  Charlestown,  United  States,  4 ARC  Center  for  Engineered  Quantum  Systems, 

School  of  Physics,  University  of  Sydney,  Sydney,  NSW,  Australia 

PURPOSE:  Free-radical-sensitive  Overhauser-enhanced  MRI  (OMRI)  is  a  promising  technique  for  imaging  the  distribution  and  dynamics  of  free 
radicals,  and  a  recently  developed  fast  high-resolution  OMRI  methodology  (1)  offers  new  perspectives  for  the  imaging  of  free  radicals  in  living  or¬ 
ganisms.  This  method  has  been  used  to  probe  BBB  breakdown  following  ischemic  stroke  in  rats  in  conjunction  with  an 
injected  stable  free  radical  (2).  Here  we  explore  if  OMRI  can  be  used  to  acquire  free  radical  images  with  sufficient  spatial 
and  temporal  resolution  to  probe  oxidative  stress  status  in  the  brain.  We  present  time-resolved  in  vivo  measurements  fol¬ 
lowing  a  single  injection  of  TEMPOL.  TEMPOL  (4-hydroxy-TEMPO),  a  small  molecule  with  a  stable  unpaired  electron 
spin,  is  detected  by  OMRI  with  very  high  sensitivity.  According  to  the  literature,  TEMPOL  does  not  cross  the  BBB  in  a 
normal  physiological  state  but  does  pass  into  brain  parenchyma  under  pathological  circumstances  associated  with  oxida¬ 
tive  stress  and  disease  (3).  As  TEMPOL  reduction  has  been  used  to  non-invasively  monitor  tissue  redox  status  in  animal 
models  (4),  we  hypothesize  that  time-resolved  OMRI  may  be  a  tool  to  elucidating  the  redox  status  of  brain  tissue  in  neuro¬ 
logic  diseases  in  which  oxidative  stress  plays  a  significant  role  such  as  reperfusion  injury  (5),  head  trauma  (6),  Alzheimer’s 
dementia  (7)  and  multiple  sclerosis  (8). 

METHODS:  A  custom  built,  low-field  MRI  scanner  equipped  with  a  biplanar  6.5  mT  electromagnet  and  biplanar  gradi¬ 
ents  was  used  in  these  experiments  (9).  3D  OMRI  was  performed  using  an  optimized  sequence  based  around  b-SSFP  as 
described  in  (1).  In  vivo  experiments  were  performed  using  a  modified  Alderman-Grant  ESR  coil  outside  of  a  high  sensi¬ 
tivity  solenoid  NMR  coil  wound  from  5/39/42  Litz  wire  that  fits  the  rat  head  (Figure  1).  Under  isoflurane  anesthesia,  the 
right  carotid  bifurcation  of  a  male  Sprague  Dawley  rat  was  dissected  and  the  pterygopalatine  artery  ligated.  Retrograde 
cannulation  (Micro-Renathane  MRE40,  Braintree  Scientific)  of  the  right  carotid  bifurcation  was  performed  through  the  external  carotid  artery.  The 
catheter,  filled  with  saline  0.9%  with  heparin  50U/ml,  was  then  tunnelized  to  the  back  of  the  animal  and  connected  to  the  injection  system.  Oxygen 
saturation,  temperature,  and  cardiac  and  respiratory  rates  were  continuously  monitored. 


start  TEMPOL  Injection  I  mL/2min.  stop 


Figure  2:  Dynamics  of  one  slice  from  3D  OMRI  dataset  acquired  from  a  rat  at  6.5  mT:  OMRI  (color)  magnitude  and  phase  images  (NA=1)  are  superimposed  on  ana¬ 
tomical  MRI  (grey),  acquired  in  the  OMRI  scanner  with  DNP  pulses  disabled  (NA=30).  Anatomical  imaging  time  was  ~5  min.  All  images,  voxel  size  l.lxl. 6x8  mm3, 
TE/TR=25/50  ms,  Matrix  128x35x11.  TEMPOL  (1.4  pl/gbw)  was  injected  into  the  right  ICA  over  2  minutes. 


Following  acquisition  of  two  reference  scans  (one  with  EPR  off  and  one  with  EPR  on),  1  ml  of  150  mM  TEMPOL  was  injected  over  2  min  with  an 
infusion  pump  (GenieTouch  Kent  Scientific).  OMRI  imaging  began  at  the  same  time  as  the  infusion,  and  continued  for  105  seconds  after  the  injec¬ 
tion  was  complete.  A  full  11  slice  OMRI  acquisition  was  acquired  every  9  seconds,  followed  by  a  6  second  delay.  This  was  repeated  16  times.  Total 
imaging  time  was  240  s. 

RESULTS:  The  internal  carotid  artery  is  responsible  of  the  vascularization  of  its  ipsilateral  hemi¬ 
sphere  via  its  terminal  branches.  The  Circle  of  Willis  enables  some  cross  vascularization  of  one 
hemisphere  from  the  contralateral  side.  Therefore,  one  would  expect  the  cerebral  distribution  of  a 
compound  injected  directly  into  the  ICA  to  reflect  the  ipsilateral  hemisphere  and  to  a  lesser  degree 
the  contralateral  hemisphere.  This  is  consistent  with  our  dynamic  in  vivo  OMRI  results  shown  in 
Figure  2  and  Figure  3,  where  marked  DNP  enhancement  is  seen  in  the  hemisphere  ipsilateral  to  the 
TEMPOL  injection.  As  the  Overhauser-enhanced  signal  has  phase  opposite  to  that  of  the  thermal 
signal,  the  phase  of  the  OMRI  image  in  Figure  2  provides  very  sensitive  contrast  even  in  cases 
where  the  radical  concentration  is  very  low  and  the  Overhauser  enhancement  may  be  small. 

DISCUSSION  &  CONCLUSION:  We  have  imaged  the  free  radical  TEMPOL  with  time-resolved 
OMRI  methods  in  vivo.  The  use  of  fast  bSSFP-OMRI  in  conjunction  with  an  exogenously  adminis¬ 
tered  free  radical  molecule  to  probe  redox  status  is  a  novel  approach.  As  a  method,  temporally  re¬ 
solved  OMRI  may  be  used  to  study  redox  status  of  brain  tissue  made  permeable  to  TEMPOL  from 
oxidative  stress.  We  note  here  that  the  long  clearance  times  seen  in  some  regions  of  the  brain  in 
Figure  2  may  indicate  TEMPOL  in  the  brain  parenchyma. 

REFERENCES:  (1)  Sarracanie  M  et  al.  MRM  2014;71(2):735-45  ;  (2)  Rosen  M  S  et  al,  Proc.  ISMRM  2014  (22)  6461  ;  (3)  Behringer  W  et  al.  J 
Cereb  Blood  Flow  Metab.  2002  22(1):  105-17  ;  (4)  Hyodo  F  et  al.  J  Pharm  Pharmacol. 2008  60(8):  1049-60  ;  (5)  Chan  PH.  Stroke.  1996  27(6):  1124-9 
;  (6)  Awasthi  D,  et  al.  Surg.  Neuro.  1997  Jun;47(6):575-81  ;  (7)  Pratico  D,  et  al  Am.  J.  Med.  2000;  109(7): 577-85  ;  (8)  (8)  Fischer  MT,  et  al  Brain. 
Oxford  Univ  Press;  2012;135(3):886-99  ;  (9)  Tsai  LL  et  al.  JMR,  2008  93(2):274-85. 
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Figure  3:  Mean  OMRI  magnitude  from  two  brain 
ROIs:  ipsilateral  (red)  and  contralateral  (blue)  to  the 
TEMPOL  injection  over  the  240  s  OMRI  acquisition. 


Figure  1:  Probe  for  OMRI 
imaging  at  6.5  mT:  NMR: 
276  kHz,  ESR:  141  MHz. 
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3D  Imaging  in  a  Portable  MRI  Scanner  using  Rotating  Spatial  Encoding  Magnetic  Fields  and  Transmit  Array  Spatial 

Encoding  (TRASE) 
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Figure  1:  (a)  Y-Z  encoding  is  done  with  the  rotating  SEM 
method,  (b)  TRASE  is  done  in  the  X  direction  using  a 
switchable  linear  Be  phase  along  X. 
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Figure  2:  The  TRASE  sequence  is  repeated  for 
every  magnet  rotation  of  the  imaging  sequence. 
The  +\-  indicates  the  sign  of  TRASE  phase  slope. 


PURPOSE:  MRI  could  find  wider  applicability  if  low-cost,  portable  systems  were  available  for  siting  in 
unconventional  locations  such  as  ICUs,  physician  offices,  ambulances,  or  rural  healthcare  sites.  The 
development  of  a  previously  described  portable  (<100kg)  MRI  scanner  for  brain  imaging1  is  advanced 
with  the  implementation  of  3rd  axis  encoding  using  the  TRASE  technique2’3. 

METHODS:  Our  2D  axial  plane  imaging  without  gradient  coils  utilized  a  rotating,  45  kg,  77  mT 
inhomogeneous  Halbach  cylinder  based  magnet1.  The  built-in  field  variation  of  the  permanent  magnet 
array  is  used  as  a  rotating  spatial  encoding  magnetic  field  (rSEM).  As  the  magnet  is  rotated  around  the 
sample,  generalized  projections  onto  the  non-linear  SEM  are  acquired  as  a  spin-echo  train. 

To  enable  imaging  in  the  3rd  direction  (along  the  axis  of  the  cylindrical  magnet),  we  use  the  TRASE 
method.  TRASE  is  a  Be  encoding  method  that  requires  2  different  switchable  Be  phase  gradients  (often 
equal  and  opposite)  along  the  encoding  direction2’3.  The  ID  TRASE  sequence  (Fig.  2)  is  a  modified 
RARE  spin-echo  train  in  which  the  slope  of  the  Be  phase  gradient  switches  for  consecutive  refocusing 
pulses  (indicated  by  the  +  or  -).  The  2  phase  gradients  correspond  to  2  k-space  origins  (-ki  and  ki),  and  the 
previous  k-space  point  is  flipped  across  the  k-space  origin  of  the  refocusing  pulse  to  traverse  k-space. 

Coils.  A  TRASE  array  was  designed  that  consists  of  two  nested  cylindrical  coils  (Fig.  3).  Coil  1  is  a  short  4- 
turn  birdcage  coil  (12  rungs,  18cm  diam.,  22cm  length)4  that  produces  a  Afield  in  Y  with  a  cosine  shape  along 
X,  5iy(x)  =  \Bfxy\cos(2n  /qx).  Coil  2  is  a  10  turn  Maxwell  coil  (22cm  diam.,  18cm  length)  that  produces  a 
Be  field  in  X  with  a  sine  shape  along  X,  Bfx(x)  =  ±\B?xy\sin(2n  /qx).  The  coils  are  tuned  to  the  Larmor 
frequency  (3.29  MHz)  and  decoupled  with  a  transformer  (-20  dB  isolation).  Figure  3b,  shows  the  relative 
magnitude  of  Bix  and  Biy  from  the  2  coils  measured  with  a  pick-up  loop.  When  the  two  coils  are  driven  with 
equal  |5^xy|  magnitude  and  spatial  frequency,  ki,  they  produce  the  desired  Be  field, 

5+ (x)  =  \Bfxy\e+i 2nkiX  (uniform  magnitude  and  linear  phase).  A  switchable  180°  phase  shifter  is  added 
in  the  Maxwell  coil  RF  path,  which  applies  a  sign  change  to  every  other  refocusing  pulse.  The  pulses  with  the 
phase  shift  thus  have  a  negative  phase  gradient;  51(x)+  =  |5^.y  |e-£  2nk^x. 

Acquisition.  The  2D  rSEM  method  and  TRASE  are  performed  simultaneously  for  3D  imaging.  The 
TRASE  array  must  rotate  with  the  magnet  because  the  Y  direction  changes  as  the  magnet  rotates.  For 
every  magnet  rotation,  the  projection  onto  the  SEM  is  acquired  as  a  TRASE-modulated  16-echo  spin- 
echo  train  (each  echo:  readout  =  196  pt,  BW  =  40  KHz,  echo  spacing  =  9.7  ms).  Data  is  acquired  at  180 
rotations  of  the  magnet  spaced  1°  apart.  Frequency-swept  WURST  pulses  are  used  to  achieve  the  same 
flip  angle  across  the  inhomogeneous  field  (BW  =  25  KHz,  [pw90,  pw180]  =  [6,3]  ms) 5.  The  birdcage  coil 
of  the  TRASE  array  is  used  as  the  receiver  coil. 

Reconstruction.  Step  one  separates  the  echo  train  data  along  X  to  obtain  16  different  YZ-datasets 
corresponding  to  each  slice.  For  each  of  the  196  points  in  the  readout  window,  the  set  of  16  points  along 
the  echo  train  populates  a  k-space  vector  (Fig.  2).  The  points  are  rearranged  and  FFT’ed  to  generate  a 
ID  projection  along  X  (Fig.  4).  Only  the  k-space  lines  formed  from  the  center  points  of  the 
echos  result  in  ID  projections  that  are  not  modulated  by  the  Bo  SEM.  This  is  repeated  for  the 
echo  trains  from  each  of  the  181  magnet  rotations.  Step  two  reconstructs  each  set  of  181 
echoes  (YZ-data)  into  a  2D  image  of  the  corresponding  slice.  The  encoding  matrix  for  each 
slice  is  calculated  using  2D  Bo  field  maps  and  then  iterative  matrix  inversion1  is  performed  to 
separately  obtain  each  2D  image. 

RESULTS:  Figure  4  shows  experimental  3D  imaging  results  using  a  7.5cm  thick  phantom  with 
1.5  cm  water-filled  compartments  spaced  3  cm  apart  in  X.  A  ID  projection  along  X  shows  the 
water-filled  compartments.  Three  slices  are  reconstructed  from  data  at  X  =  -3cm,  Ocm,  and  3cm 
using  the  corresponding  2D  Bo  field  maps. 

DISCUSSION:  The  non-bijectivity  of  the  SEMs  results  in  image  aliasing  across  the  center.  This 
is  usually  resolved  using  a  surface  coil  array6,  but  as  a  first  step,  the  birdcage  coil  of  the  array 
was  used  for  a  receiver  coil  and  only  the  reduced  FOV  image  is  shown.  The  blurring  in  the 

image  slice  is  due  to  a  combination  of  the  spatially  varying  resolution  of  the  SEMs  and  Figure  4:  Image  of  a  3D  phantom  consisting  of  three  water-filled 

systematic  calibration  errors  in  our  model-based  reconstruction.  compartments  spaced  3  cm  apart  in  X.  A  ID  projection  along  X 

CONCLUSION:  We  demonstrate  proof-of-concept  3D  images  in  an  inhomogeneous  45  kg  for  1  an9le,'s  shown  ('ower  left^  3  image  slices  were 

magnet  without  the  use  of  the  gradient  coils.  TRASE  encoding  along  X  is  ideal  because  like  our  reconstructed  usmg  the  appropnate  2D  field  maps. 

2D  method  it  uses  spin-echo  trains,  and  it  prevents  the  need  for  a  Bo  gradient  coil  and  GPA,  which  would  add  weight,  cost,  and  power  requirements  to  the  portable, 
low-cost  system.  Future  work  will  extend  this  3D  imaging  method  using  larger  diameter  magnets  and  coils  with  the  goal  of  portable  low-cost  MRI  brain  imaging. 
REFERENCES:  (1)  Cooley  C,  MRM  2014.  (2)  Sharp  JC,  MRM  2010.  (3)  Sharp  JC,  NMR  in  Biomed.  2013.  (4)  Borsboom  H,  Magn.  Reson.  Mater.  Phy.  1997.  (5) 
Casabianca  LB,  JMR  2014,  (6)  Hennig  J,  Magn.  Reson.  Mater.  Phy,  2008. 
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Figure  3:  (a)  The  TRASE  array  is  a  nested  birdcage 
and  Maxwell  coil,  (b)  The  measured  magnitude  of  Bu 
from  the  Maxwell  coil  and  Biy  from  the  birdcage  coil. 
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High  Performance  Probe  for  in  vivo  Overhauser  MRI 
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Purpose:  Overhauser-enhanced  MRI  (OMRI)  is  an  electron-proton  double  resonance  imaging  tech¬ 
nique  of  much  interest  due  to  its  ability  to  detect  the  concentration  and  distribution  of  free  radicals. 
Tracking  of  exogenous  free  radicals  with  OMRI  in  vivo  has  enabled  the  development  of  oxymetry 
probes  [1]  and  the  imaging  of  redox  reactions  [2].  The  large  gyromagnetic  ratio  of  electrons  (28 
GHz/T)  demands  that  in  vivo  OMRI  is  performed  at  very  low  magnetic  fields  (~10  mT)  in  order  to 
minimize  RF  heating  and  penetration  depth  issues.  Operation  at  low  magnetic  field  causes  a  drastic 
reduction  in  NMR  sensitivity  despite  the  signal  enhancement  that  comes  from  the  Overhauser  effect, 
and  emphasizes  the  need  for  high  S/N  probes.  OMRI  probe  design  is  still  relatively  unexplored, 
despite  its  importance,  and  presents  challenges  unique  to  the  frequencies  of  operation  (fH  =  276  kHz 
and  fe  =  140.8  MHz  in  our  experiments  at  6.5  mT).  Here,  we  report  the  development  of  a  high  per¬ 
formance  OMRI  probe  built  to  image  free  radical  probes  of  the  blood  brain  barrier  following  ischem¬ 
ic  stroke  in  a  rat  model  [3]. 

Methods:  Our  OMRI  probe,  consisting  of  an  NMR  solenoid  inside  a  modified  Alderman-Grant 
Resonator  (Figure  1),  was  designed  for  use  with  a  custom  built,  very  low-field  MRI  scanner  operat¬ 
ing  at  6.5  mT  (fH  =  276  kHz)  [4],  Images  were  acquired  using  a  recently  developed,  fast,  high- 
resolution  b-SSFP  based  OMRI  methodology  [5].  NMR  probe  design  in  the  low-field  regime  is 
fundamentally  different  to  that  at  conventional  MRI  fields  as  thermal  noise  due  to  the  intrinsic  re¬ 
sistance  of  the  pickup  coil  dominates  over  sample  noise.  This  leads  to  a  compromise  where  S/N 
improvements  come  at  the  expense  of  imaging  bandwidth  S/N  ~y[Q  ~1  /s/BW.  An  85  turn  solenoid 
was  wound,  using  low  AC  resistance  5/39/42  litz  wire,  on  a  3D  printed  polycarbonate  former.  This 
high  filling  factor  coil  has  a  bandwidth  of  3  kHz. 


Figure  2:  COMSOL  Multiphysics  simulations  (a)  Our  modified  Alderman-Grant  Resonator,  (b) 
Plot  of  |E|  showing  that  the  electric  field  is  strongly  suppressed  inside  the  resonator,  (c)  Bie  shows 
high  homogeneity  in  the  imaging  volume  -  less  than  10%  variation  across  the  imaging  region. 
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Figure  1:  NMR/ESR  rat  head  probe  for  operation  at  6.5  mT. 
The  ESR  resonator  is  tuned  to  140.8  MHz.  The  litz  wire 
NMR  solenoid  coil  (276  kHz)  resonator  board  is  not  shown. 


A  modified  Alderman-Grant  ESR  resonator  was  built  using 
copper  foil  on  Pyrex  tubing.  All  metal  placed  in  close  prox¬ 
imity  to  the  NMR  solenoid  strongly  couples,  reducing  the 
NMR  sensitivity.  We  therefore  minimized  amount  of  copper 
in  the  ESR  resonator.  Windows  were  removed  from  the  pan¬ 
els  on  the  sides  of  an  Alderman-Grant  resonator,  a  region  of 
low  current  flow  [6],  in  an  attempt  to  reduce  coupling  whilst 
maintaining  Bi  homogeneity.  Shielding  at  the  ends  of  the 
resonator  prevents  high  electric  fields  at  the  capacitors  pene¬ 
trating  the  imaging  volume,  important  because  PRF-absorbed  ~ 
E2.  Slits  in  the  shielding  prevent  the  formation  of  closed 
loops  that  couple  to  the  solenoid.  Figure  2  demonstrates  the 
high  Bi  homogeneity  and  strong  E  suppression  in  our  ESR 
resonator. 


TEMPOL  (4-hydroxy-TEMPO)  a  stable  radical,  is  detected  by  OMRI  with  high  sensitivity,  and 
because  of  its  small  size  may  facilitate  imaging  blood  brain  barrier  leakage  in  cases  of  oxidative 
stress  [3].  Simulations  were  validated  in  2mM  TEMPOL  solutions  (not  shown)  using  simple 
spectroscopic  measurements  as  well  as  fast  imaging  strategies  [5].  The  ESR  resonator  was  tuned 
to  141  MHz,  the  lowest  frequency  of  the  TEMPOL  triplet  state  at  6.5  mT,  to  minimize  RF  heating 
during  ESR  irradiation. 

Results:  Previously  we  reported  the  development  of  an  OMRI  probe  for  a  rat  head  model  utilizing 
an  NMR  solenoid  and  ESR  surface  coil  [3].  Testing  of  our  new  OMRI  probe  shows  that  it  has  3x 
the  S/N  of  the  probe  reported  in  [3]  and  rectifies  problems  with  Bie  homogeneity,  yielding  homo¬ 
geneous  enhancement  of  -6.7  in  2mM  TEMPO  when  10  W  of  RF  power  is  applied. 

In  vivo  OMRI  signal  enhancement  is  clearly  visible  in  the  rat  brain  after  TEMPOL  injection,  as 
shown  in  Figure  3.  As  the  Overhauser-enhanced  signal  has  a  phase  opposite  to  that  of  the  thermal 
signal,  the  phase  image  in  Figure  3  provides  sensitive  contrast  in  regions  of  low  TEMPOL  con¬ 
centration. 

Discussion  and  Conclusion:  We  have  implemented  a  high  performance  probe  for  high  temporal 
and  spatial  resolution  OMRI  in  a  rat  brain.  The  S/N  of  this  probe  may  be  further  improved,  whilst 
maintaining  imaging  bandwidth,  by  using  an  active  feedback  circuit  [7].  Higher  S/N  could  also 
be  realized  through  supercooling  of  the  NMR  solenoid  or  by  using  a  free  radical  with  a  longer  T  k. 
such  as  triphenylmethyl  [1],  This  probe  may  allow  the  in  vivo  detection  of  rapid  redox  changes  in 
pathologic  tissues,  specifically  in  the  context  of  brain  trauma  or  stroke. 

References:  [1]  K.  Golman  et  al.,  J.  Mag.  Res.  Im„  12,  929-938  (2000).  [2]  Utsumi  et  al.,  PNAS,  103,  1463-1468  (2006).  [3]  M.  Rosen  et  al„  ISMRM,  6461  (2014). 
[4]  Tsai  et  al.,  J.  Mag.  Res.,  193,  274-285  (2008).  [5]  M.  Sarracanie  et  al.,  Magn  Res.  Med.,  71,  735-745  (2014).  [6]  Mispelter  et  al.  ‘NMR  Probeheads:  For  Biophysical 
and  Biomedical  Exeperiments’,  Imperial  College  (2006).  [7]  Baudin  et  al.,  J.  of  Phys.  Conf.  Series,  294,  012009  (2011). 
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Figure  3:  OMRI  images  acquired  from  a  rat  at  6.5  mT  follow¬ 
ing  injection  of  lmL  of  150  mM  TEMPOL.  Five  slices  from 
an  1 1  slice  data  set  shown.  OMRI  (NA=1)  imaging  time  was  9 
seconds.  Anatomical  MRI  (NA=30)  was  acquired  in  the  OMRI 
scanner  with  ESR  power  disabled.  All  images,  voxel  size:  1.1 
x  1.6  x  8  mm3,  Matrix:  128  x  35  x  11 
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Purpose:  MRI  at  low  magnetic  field  (<10  mT)  without  cryogenic  or  hyperpolarization  techniques  presents  unique  engineering  challenges.  Imaging 
coils  must  maximize  coverage  over  the  volume  of  interest  while  minimizing  losses  in  a  regime  unusual  in  contemporary  MRI — where  Johnson  noise 
dominates  the  noise  floor.  Our  previous  work  validated  our  scanner  hardware  and  pulse  sequences  MRI  at  6.5  mT.  The  purpose  of  the  present  work  is 
to  construct  a  high-performance  coil  for  low-field  imaging  of  the  human  brain  in  vivo. 

Methods:  Our  previous  coils  focused  on  implementing  array  coil  technology  common  at  clinical  field  strengths,  and  we  constructed  an  8  channel 
receive-only  coil  for  our  6.5  mT  low-field  imager  at  276  kHz.  While  our  results  were  promising1’2,  a  lack  of  low  noise  high-impedance  pre-amplifiers 
at  our  frequency  regime  prevented  us  from  achieving  sufficient  coil  decoupling,  resulting  in  images  with  poor  SNR.  We  present  here  a  single  channel 
coil  with  greatly  improved  SNR  for  these  applications. 

We  designed  a  spiral  volume  coil  design  for  combined  Tx  and  Rx  and  optimized  for  human  head  imaging  in  the  trans¬ 
verse  BO  field  of  our  low  field  scanner.  This  places  the  receive  element  close  to  the  head  while  still  covering  the  entire 
volume  uniformly.  Additionally,  the  symmetry  of  a  spherical  spiral  results  in  a  homogeneous  magnetic  field  over  its 
volume,  minimizing  B1  inhomogeneity3’4.  Losses  in  low  field  imaging  are  mainly  due  to  resistive  losses  in  the  coil 
(the  so-called  Johnson  noise  regime).  Rather  than  lower  this  resistance  by  cooling,  we  instead  use  of  multistranded 
Litz  wire.  When  compared  to  a  solid  copper  wire  of  the  same  gauge,  a  loop  made  of  Litz  wire  will  have  the  same  in¬ 
ductance  but  a  fraction  of  the  AC  resistanc,  resulting  in  a  lower  noise  floor  and  higher  image  SNR. 

The  optimal  number  of  turns  and  appropriate  wire  gauge  were  determined  through  simulation.  The  helmet  shape  was 
designed  to  fit  closely  to  an  anatomically  accurate  human  head  model.  The  spiral  path  was  generated  in  a  custom 
MATLAB  (Mathworks  Inc.)  script  and  imported  into  BiotSavart  (Ripplon  Software  Inc.)  for  an  estimate  of  induct- 
Figure  1:  Form  fitting  helmet  ance.  Resistance  was  computed  from  the  estimated  wire  length  and  the  characteristic  impedance  per  length  of  wire  for 
with  30-tum  spiral  design.  that  gauge.  Coil  Q  was  estimated  as  coL/R.  Bandwidth  was 
estimated  from  the  Q  value.  We  iterated  through  several 
designs  until  we  achieved  our  target  Q~40  and  BW=6kHz. 

Results:  From  simulations,  we  built  a  30-turn  spiral  wound  with  Type  1  40/38  Litz  wire. 

Our  aligned  tum-to-tum  spacing  was  5.6  mm.  The  final  design  was  3D  printed  using 
fused  deposition  modeling  technology  with  polycarbonate  in  a  Fortus  360  me  printer 
(Stratasys,  Eden  Prairie,  MN,  USA),  (Figure  1). 

Images  were  acquired  at  6.5  mT  in  a  head  shaped  phantom  using  3D  b-SSFP  with  50% 
undersampling.  Resulting  voxel  size  was  3.9x3.5x16.7  mm3  (NA=110)  (Figure  2a).  A 
maximum  in-plane  SNR  of  130  was  calculated.  An  image  of  a  similar  slice  obtained  with 
the  8-channel  array  coil  is  shown  for  comparison  (Figure  2b,  combined  image  using  RSS, 
voxel  size  =  4.4x4.2x22  mm3,  NA  =  80,  max  SNR  =  16.5). 

19  slices  of  structured  resolution  phantom  imaged  with  optimized  spiral  volume  coil  are 
shown  in  Figure  3.  This  four-quadrant  phantom  is  filled  will  liquid-filled  spheres  of  various  diameters.  Voxel  size  in  this  scan  is  3.5x2. 5x5  mm3  with 
NA=40  and  maximum  SNR  of  24.8.  Total  acquisition  time  is  10  minutes. 


Figure  2:  Equivalent  slices  in  a  head  phantom  acquired  using 
the  single  channel  spiral  coil  (A)  and  the  8 -channel  array  (B  ) 
at  6.5  mT.  Images  are  scaled  by  their  maximum  intensity. 


Figure  3:  Axial  view  of  a  multi-compartment  head-shaped  phantom  acquired  with  the  Litz  spiral  volume  coil  at  6.5  mT  in  10  minutes. 


Discussion:  A  30-turn  single  channel  volume  spiral  coil  designed  for  low  field  imaging  significantly  outperformed  our  8-channel  array  coil.  By  min¬ 
imizing  losses  in  the  coil,  maximizing  filling  factor  and  eliminating  coupling  issues,  we  were  able  to  obtain  significantly  higher  SNR. 

Conclusions:  A  high  filling  factor  human-head  Tx/Rx  spiral  volume  imaging  coil  was  successfully  designed  and  constructed  for  6.5  mT  (276  kHz). 

Financial  Disclosure:  This  research  was  supported  by  the  Department  of  Defense,  Defense  Medical  Research  and  Development  Program,  Applied 
Research  and  Advanced  Technology  Development  Award  W81XWH- 11-2-0076  (DM09094). 
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Transmit  Array  Spatial  Encoding  (TRASE)  with  broadband  WURST  pulses  for  robust  spatial  encoding  in  inhomogeneous  BO 

fields 
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TARGET  AUDIENCE:  MRI  engineers  working  with  inhomogeneous  B0  fields 
PURPOSE:  Transmit  array  spatial  encoding  (TRASE)1  uses  RF  coils  with  linearly-varying 
Bf  phase  to  replace  B0  gradient  phase  encoding  along  one  more  directions.  To  traverse  k- 
space,  the  sign  of  the  Bf  phase  is  flipped  for  alternating  refocusing  pulses  in  a  spin  echo 
train,  imparting  an  additional  Ak  phase  modulation  to  each  successive  echo  in  the  train. 
TRASE  provides  an  alternative  for  performing  high  resolution  spatial  encoding  in 
applications  where  conventional  B0  gradients  are  not  feasible.  However,  the  TRASE  point 
spread  function  degrades  rapidly  for  refocusing  pulse  angles  outside  the  range  150-195° 
due  to  accumulation  of  flip  angle  and  phase  errors  along  the  echo  train2.  For  the  hard  pulses 
(as  shown  in  TRASE  results  to  date),  the  refocusing  angle  falls  below  150°  when  AB0  >  40% 
of  the  pulse  Rabi  frequency.  This  constrains  the  use  of  TRASE  in  applications  where  peak 
RF  power  is  limited  or  AB0  is  large.  We  propose  a  method  for  overcoming  this  problem  using 
broadband  “WURST”  RF  pulses  along  with  a  simple  deconvolution  method  to  render  the 
acquired  echoes  compatible  with  TRASE.  METHODS:  To  demonstrate  ID  TRASE  in  the 
presence  of  strong  off-resonance,  we  acquire  projections  in  an  inhomogeneous  Halbach 
array  magnet  (±20  KHz  over  FOV)  originally  developed  as  a  lightweight,  permanent  magnet- 
based  portable  MRI  prototype3  with  Larmor  frequency  3.285  MHz  (Fig.  1).  WURST-40 
pulses4  with  phase  cycling5  are  used  for  spin  excitation  and  refocusing.  Pulses  with  linear 
frequency  sweep  impart  a  quadratic  phase  across  the  spin  bandwidth,  but  if  the  pulse  sweep 
rate  satisfies  the  condition  Rexc  =  VzRret  and  the  AB0  field  is  constant,  then  the  refocusing 
pulse  removes  the  phase  applied  by  the  excitation  pulse  and  all  isochromats  refocus  at  7"E\ 
Further  readouts  in  the  echo  train  alternate  between  conventional  echoes  and  “spectral” 
echoes5  that  resemble  the  frequency  distribution  of  spins  in  the  AB0  field.  This  additional  RF 
phase  modulation  must  be  removed  in  order  to  make  the  spectral  echoes  compatible  with 
TRASE.  The  spectral  echo  signal  can  be  approximated  as  a  convolution  of  the  conventional 
spin  echo  signal  with  a  chirp  function,  exp(i27r  t 2),  where  depends  on  the  sweep  rate. 
DECONVOLUTION  is  performed  by  either  (a.)  convolving  the  spectral  echo  with  a  scaled, 
conjugated  version  of  the  chirp  function,  or  equivalently  (b.)  in  the  Fourier  transform  (FT) 
domain  by  dividing  the  FT  of  the  echo  (ID  projection)  by  the  FT  of  the  chirp  function7.  Since 
the  FT’ed  chirp  kernel  is  a  pure  phase  term,  method  (b.)  requires  simply  taking  the 
difference  between  the  phase  of  the  first  and  second  projections  in  the  echo  train  (Fig.  2). 
This  phase  difference  is  used  to  correct  all  subsequent  spectral  echoes.  The  FT  of  the 
corrected  spectral  echoes  yields  echo  magnitudes  and  phases  that  correspond  closely  to  the 
conventional  echoes  in  the  train.  EXPERIMENT:  Tube  phantoms  are  placed  in  the 
inhomogeneous  B0  field  and  a  train  of  129  echoes  are  acquired  each  with  196  readouts 
points,  40  KHz  BW,  and  echo  spacing  9.4ms.  WURST-40  excitation  and  refocusing  pulse 
durations  are  [  exc,  ref]  =  [6,  3]  ms  and  sweep  rates  are  [Rexc,  Rret ]  =  [4.15,  8.3]  MHz/s  over  a 
25  KHz  bandwidth.  A  multi-turn  Maxwell  coil  ( B1x+ )  is  paired  with  a  concentric  multi-turn 
birdcage  coil8  ( B1y+ )  to  generate  sine  and  cosine-shaped  Bf  amplitudes,  respectively,  along 
the  encoding  axis,  x.  When  both  coils  are  driven  simultaneously,  the  desired  linear  phase 
slope  is  generated:  B^tot(x)  =  \B^xy\cos(27rkx)  ± i\B^xy\sin(2wkx)  =  \B+Xy\exp (±i2Trkx) 

The  TRASE  Bf{x)  phase  is  switched  using  an  RF  phase  shifter  to  play  the  pattern: 
{90+,180+, echo,  180+,ec/7O,180-,ec/7O,180+,ec/70, 180-...}.  The  first  two  echoes  have  no 
TRASE  modulation  and  can  therefore  be  used  to  calculate  the  phase  difference  required  for 
the  spectral  echo  correction.  To  create  a  ID  projection  along  x,  the  point  from  the  center  of 
each  readout  is  selected  to  avoid  modulation  imparted  by  the  AB0  field  (Fig.  3).  The  points  are 
reordered  and  FT’ed  to  create  a  ID  projection  of  the  object.  RESULTS:  Fourier  domain  phase 
subtraction  faithfully  restores  spectral  echoes  so  that  their  real  and  imaginary  components 
align  closely  with  those  of  conventional  echoes.  The  TRASE-WURST  acquisition  successfully 
generates  projections  of  phantoms  with  bandwidths  >10  KHz  in  the  inhomogeneous  B0  field. 
By  contrast,  recognizable  projections  could  not  be  obtained  using  hard  pulse  excitations  since 
available  RF  power  constrained  the  minimum  pulse  duration  ([75,150]  ps  for  [90°,  180°]  pulses) 
and  thus  the  bandwidth  of  the  pulses.  CONCLUSION:  We  demonstrate  robust  spatial 
encoding  using  WURST  pulses  for  TRASE  encoding  in  an  inhomogeneous  field.  We  further 
show  a  simple  phase  correction  that  allows  WURST  spectral  echoes  to  be  used  for  TRASE 
encoding.  In  future  work,  TRASE-WURST  slice  encoding  will  be  combined  with  in-plane  B0 
gradient  encoding  to  generate  3D  images.  Potential  applications  include  silent  imaging,  oil  well 
prospecting,  and  portable  low-field  scanners.  REFERENCES:  [1]  Sharp  JC,  NMR  Bio.  2013. 
[2]  Sharp  JC,  ISMRM  2010,  p.  1469.  [3]  Cooley  C,  MRM  2014.  [4]  Kupce  E,  JMR  1995.  [5] 
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Fig.  1.  Halbach  magnet  AB0  field  plotted  along  three 
planes.  TRASE  encoding  is  implemented  along  x direction. 
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Fig.  2.  WURST  pulse  sequence  diagram  with  four  echoes 
acquired  on  a  water  bottle  phantom.  The  bracketed  block  is 
repeated  as  needed.  Conventional  echoes  alternate  with 
“spectral”  echoes  modulated  by  the  quadratic  phase  of  the 
frequency  swept  RF  pulses.  This  effect  is  compensated  by 
removing  the  quadratic  phase  from  the  projections 
corresponding  to  each  spectral  echo.  After  phase  correction 
the  echo  real  and  imaginary  parts  overlay  for  all  echoes  in 
the  train.  Phase  cycling  scheme  is  as  follows:  </>  =[90,  270, 
180,  0]°,  <t>re=\ 0,  0,  90,  90]°,  </)^f=[90,  90,  0,0]°,  =[0, 
180,  90,  270]°.  The  3rd  and  4th  echoes  in  the  block  are 
multiplied  by  it  to  compensate  for  phase  cycling. 
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Fig.  3.  The  first  8  corrected  echoes  (a.)  from  a  128-echo 
train  shown  with  (red)  and  without  TRASE  phase 
modulation.  Orange  and  green  dots  at  each  echo  time 
denote  positive  and  negative  k-space  points,  respectively, 
which  are  reordered  into  ID  k-space  (b.).  The  Fourier 
transform  yields  a  128-point  ID  projection  (c.)  of  the 


water-filled  tube  phantom  (d.). 
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Two-Dimensional  Imaging  in  a  Lightweight  Portable 
MRI  Scanner  without  Gradient  Coils 
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Purpose:  As  the  premiere  modality  for  brain  imaging,  MRI 
could  find  wider  applicability  if  lightweight,  portable  systems 
were  available  for  siting  in  unconventional  locations  such  as 
intensive  care  units,  physician  offices,  surgical  suites,  ambu¬ 
lances,  emergency  rooms,  sports  facilities,  or  rural  healthcare 
sites. 

Methods:  We  construct  and  validate  a  truly  portable  (<100 
kg)  and  silent  proof-of-concept  MRI  scanner  which  replaces 
conventional  gradient  encoding  with  a  rotating  lightweight 
cryogen-free,  low-field  magnet.  When  rotated  about  the 
object,  the  inhomogeneous  field  pattern  is  used  as  a  rotating 
spatial  encoding  magnetic  field  (rSEM)  to  create  generalized 
projections  which  encode  the  iteratively  reconstructed  two- 
dimensional  (2D)  image.  Multiple  receive  channels  are  used  to 
disambiguate  the  nonbijective  encoding  field. 

Results:  The  system  is  validated  with  experimental  images  of 
2D  test  phantoms.  Similar  to  other  nonlinear  field  encoding 
schemes,  the  spatial  resolution  is  position  dependent  with 
blurring  in  the  center,  but  is  shown  to  be  likely  sufficient  for 
many  medical  applications. 

Conclusion:  The  presented  MRI  scanner  demonstrates  the 
potential  for  portability  by  simultaneously  relaxing  the  magnet 
homogeneity  criteria  and  eliminating  the  gradient  coil.  This  new 
architecture  and  encoding  scheme  shows  convincing  proof  of 
concept  images  that  are  expected  to  be  further  improved  with 
refinement  of  the  calibration  and  methodology.  Magn  Reson 
Med  000:000-000,  2014.  ©  2014  Wiley  Periodicals,  Inc. 
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INTRODUCTION 

Specialized,  portable  MRI  systems  have  the  potential  to 
make  MR  neuroimaging  possible  at  sites  where  it  is  cur¬ 
rently  unavailable  and  enable  immediate,  “point-of-care” 
detection  and  diagnosis  of  acute  intracranial  pathology 
which  can  be  critical  in  patient  management.  For  exam¬ 
ple,  the  characterization  of  acute  posttraumatic  space 
occupying  brain  hemorrhage  is  a  time-sensitive  emer¬ 
gency  for  which  simple  clinical  assessment  and  even 
urgent  computed  tomography  scanning  may  be  insuffi¬ 
cient.  While  conventional  MR  scanners  are  capable  of 
making  this  diagnosis,  they  are  not  available  in  remote 
locations.  In  Intensive  Care  Units  scanners  are  generally 
nearby,  but  are  difficult  to  use  because  of  the  dangers 
associated  with  transporting  critical  care  patients.  A 
portable  bed-side  scanner  could  offer  major  benefits  in 
such  situations.  Portable,  low-cost  scanners  are  compel¬ 
ling  for  applications  where  power,  siting,  and  cost  con¬ 
straints  prohibit  conventional  scanners.  Examples 
include  clinics  in  rural  or  underdeveloped  areas,  mili¬ 
tary  field  hospitals,  sports  arenas,  and  ambulances. 
Finally,  analogous  to  the  current  use  of  ultra-sound,  a 
low-cost  and  easy-to-implement  scanner  could  find  uses 
in  neurology,  neurosurgery,  or  neuro-oncology  examina¬ 
tion  rooms  for  routine  disease  monitoring  (e.g.,  monitor¬ 
ing  ventricle  size  after  stent  placement).  The 
development  of  a  portable  scanner  relies  on  the  co¬ 
design  of  a  new  image  encoding  methods  and  simplified 
hardware.  This  approach  is  detailed  in  the  present  work. 

Traditional  Fourier  MR  imaging  methods  rely  on 
homogeneous  static  polarizing  fields  ( B0 )  and  high 
strength  linear  spatial  encoding  magnetic  fields  (SEMs) 
produced  by  magnetic  gradient  coils.  Conventional  scan¬ 
ners  use  high  cost  superconducting  wire,  liquid  cryogen 
cooling  systems,  and  high  power  supplies  and  electron¬ 
ics.  These  aspects  make  it  difficult  to  simply  scale-down 
conventional  MRI  scanners  to  portable,  low  cost  devices. 
Recently  in  Kose  and  Haishi  (1)  and  Kimura  et  al  (2), 
high  resolution  imaging  has  been  achieved  with  table- 
top  and  small  bore  permanent  magnet  systems  with  long 
acquisition  times  (1),  including  a  mobile  MRI  system 
developed  for  outdoor  imaging  of  small  tree  branches  (2), 
but  these  scanners  lack  a  bore  size  suitable  for  brain 
imaging  and  the  long  acquisition  times  are  not  conducive 
to  imaging  in  triage  settings.  Other  approaches  that 
scale-down  the  size  of  conventional  systems  for  intra¬ 
operative  MRI  show  promise  (3).  However,  while  these 
systems  are  relatively  easy  to  retrofit  in  operating  room, 
they  are  not  truly  portable. 
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In  the  present  work,  we  use  a  novel  image  encoding 
method  based  on  rotating  spatial  encoding  magnetic 
fields  (rSEMs)  to  create  a  portable  scanner.  We  replace 
the  B0  magnet  and  linear  gradient  coils  with  a  rotating 
permanent  magnet  featuring  an  inhomogeneous  field  pat¬ 
tern  used  for  spatial  encoding.  In  this  scheme,  the  inho¬ 
mogeneity  in  the  B0  field  serves  as  a  spatial  encoding 
magnetic  field  (SEM),  and  is  requirement  for  image 
encoding  not  a  nuisance.  Loosening  the  homogeneity 
constraint  of  some  conventional  magnet  designs  leads  to 
a  reduction  in  the  minimum  required  magnet  material, 
and  allows  for  more  sparse/lightweight  designs  (45  kg  in 
our  prototype).  Additionally,  the  rotation  of  the  magnet’s 
inhomogeneous  field  pattern  replaces  the  function  of 
heavy  switchable  gradient  coils  with  significant  power 
requirements. 

Several  NMR  devices  for  niche  applications  have 
explored  relaxing  the  magnet  homogeneity  constraint,  as 
well  as  reducing  the  reliance  on  traditional  Fourier 
image  encoding.  The  oil  well-logging  industry  was  the 
first  to  explore  the  idea  of  mobile  NMR  using  “external 
sample”  or  “inside-out”  NMR  tools  for  measuring  fluid 
in  rock  formations  down-hole  (4).  This  work  was  ini¬ 
tially  done  with  electromagnets  or  in  the  earth’s  field, 
but  the  advent  of  rare-earth  magnets  with  high  energy 
products  such  as  SmCo  and  NdFeB  (5),  has  allowed 
more  effective  borehole  NMR  tools  to  be  developed  (6). 
Some  portable  single-sided  NMR  devices  (7,8)  exploit 
inhomogeneous  magnetic  fields  from  permanent  magnets 
for  ID  spatial  encoding.  In  these  systems  a  rare-earth 
magnet  array  is  placed  against  the  object  such  that  the 
field  falls  off  roughly  linearly  with  depth.  Broadband 
excitation  and  spin-echo  refocusing  are  used  to  obtain  a 
one-dimensional  (ID)  depth  profile  of  the  water  content 
in  the  object  (9-11).  Thus,  these  systems  use  the  inhomo¬ 
geneity  of  the  small  magnet  to  spatially  encode  the  depth 
of  the  water;  a  principle  that  we  will  exploit  in  a  more 
complete  way. 

Previously,  Cho  et  al  implemented  a  mechanically 
rotating  DC  gradient  field  in  conjunction  with  a  conven¬ 
tional  MRI  scanner  with  the  motivation  of  silent  imaging 

(12) .  In  that  case,  the  rotating  electromagnet  produced  a 
linear  gradient  field  so  traditional  projection  reconstruc¬ 
tion  methods  could  be  used.  In  the  presently  described 
portable  scanner,  the  dominant  SEM  field  term  is  quad- 
rupolar,  which  requires  specialized  acquisition  and 
reconstruction  techniques.  Spatial  encoding  with  similar 
nonlinear  fields  created  with  electromagnets  has  recently 
drawn  attention  as  a  way  to  achieve  focused  high  imag¬ 
ing  resolutions,  reduced  peripheral  nerve  stimulation 

(13) ,  and  improved  parallel  imaging  performance  (14).  In 
our  scanner,  the  approximately  quadrupolar  SEM  fields 
are  physically  rotated  around  the  object  along  with  the 
B0  field,  and  stationary  RF  coils  are  used  to  acquire  gen¬ 
eralized  projections  of  the  object  in  spin-echo  train  form. 

In  this  manuscript,  we  describe  the  design,  construc¬ 
tion,  and  testing  of  a  portable  2D  MRI  scanner.  We  show 
that  the  encoding  scheme  we  introduce  can  achieve  a 
resolution  of  a  few  millimeters  in  phantom  images. 
While  full  3D  encoding  is  not  demonstrated,  the  system 
is  compatible  with  RF  encoding  schemes,  such  as  the 
transmit  array  spatial  encoding  (TRASE)  method  (15,16), 


capable  of  adequately  encoding  the  third  dimension 
(along  the  axis  of  the  cylindrical  magnet).  The  magnet 
design  and  initial  encoding  attempts  were  previously 
reported  in  abstract  form  (17-19). 

METHODS 

Magnet  and  Field  Mapping 

The  described  rSEM  encoding  method  is  valid  for  arbi¬ 
trary  encoding  field  shapes,  although  the  shape  will 
affect  the  spatially  variable  resolution  of  the  images.  A 
sparse  dipolar  Halbach  cylinder  design  similar  to  the 
“NMR  Mandhala”  (20,21)  was  chosen  to  produce  the 
rotating  B0  field  presented  here.  This  arrangement  of  per¬ 
manent  magnets  produces  an  approximately  uniform 
field  directed  transverse  to  the  axis  of  the  cylinder  (22). 
The  major  design  criteria  for  our  Halbach  magnet  were: 
(i)  maximum  average  field  for  highest  signal  to  noise 
ratio  (SNR),  (ii)  sufficient  field  variation  for  spatial 
encoding  while  maintaining  reasonable  measurement 
and  excitation  bandwidths,  (iii)  minimum  volume  of  per¬ 
manent  magnet  material  to  keep  cost  and  weight  down, 

(iv)  use  of  stock  rare-earth  magnet  material  shapes,  and 

(v)  minimum  size  to  fit  the  head  (to  maximize  B0).  Note, 
the  design  was  not  focused  on  the  specific  spatial  encod¬ 
ing  field  shape,  and  the  resulting  pattern  in  the  con¬ 
structed  magnet  was  accepted  as  the  SEM  for  the 
presented  scanner.  The  Halbach  cylinder  design  consists 
of  a  36-cm-diameter  array  of  20  rungs  comprising  N42 
grade  NdFeB  magnets  that  are  each  1  x  1  x  14"  (magne¬ 
tized  through  the  1"  thickness).  Two  additional  Halbach 
rings  made  up  of  twenty  1"  NdFeB  cubes  are  added  to 
the  ends  of  the  cylinder  to  reduce  field  fall-off  caused  by 
the  relatively  short  length  of  cylinder.  Figure  1  shows  a 
3D  drawing  of  the  magnet,  the  simulated  field,  and  con¬ 
structed  magnet. 

The  predicted  field  pattern  (Fig.  la,b),  as  well  as  the 
forces  between  the  magnet  rungs  were  simulated  using 
COMSOL  Multiphysics  (Stockholm,  Sweden).  This  cal¬ 
culation  estimated  an  internal  force  of  178  N,  which  is 
adequately  handled  by  the  fiberglass  and  ABS  frame 
designed  to  hold  the  NdFeB  magnet  array.  The  magnet 
rungs  consist  of  NdFeB  magnets  stacked  inside  square 
fiberglass  tubes  (McMaster-Carr,  Elmhurst,  IL,  USA), 
which  are  fixed  by  five  water  jet  cut  3/8"ABS  rings  (Fig. 
Id).  Each  14"  long  magnet  rung  is  comprised  of  four 
individual  bar  magnets  (Applied  Magnet,  Plano,  TX, 
USA)  which  were  bonded  together  (three  4"  bars  and  a 
2"  bar).  The  ABS/fiberglass  frame  was  assembled  before 
NdFeB  magnet  handling,  and  then  the  magnet  rungs 
were  populated  one  at  a  time.  Because  the  four  magnets 
comprising  each  rung  repel  each  other  as  they  are 
inserted,  a  magnet  loading  and  pushing  jig  was  neces¬ 
sary  to  force  the  magnets  together  while  the  magnet 
bonding  adhesive  cured  (Loctite  p/n  331  and  7387, 
Diisseldorf,  Germany).  The  jig  was  a  simple  threaded 
rod  mounted  to  the  magnet  assembly  frame  above  the 
opening  of  the  fiberglass  tube. 

The  constructed  magnet  weighs  45  kg  and  has  a  77.3 
mT  average  field  in  the  16-cm  field  of  view  (FOV)  center 
plane,  corresponding  to  a  3.29  MHz  proton  Larmor  fre¬ 
quency.  The  cylindrical  magnet  sits  on  aluminum  rollers 
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FIG.  1.  The  magnet  array  consists  of  twenty  1"  x  1"  x  14"  NdFeB  magnets  oriented  in  the  k  =  2  Halbach  mode.  Additional  Halbach 
rings  made  of  1"  x  1"  x  1"  magnets  were  added  at  the  ends  to  reduce  field  fall  off  along  the  cylindrical  axis.  a,b:  Simulation  of  the 
magnetic  field  in  two  planes.  The  field  is  oriented  transverse  to  the  cylinder  axis  (z-direction).  c:  Schematic  of  NdFeB  magnets  compos¬ 
ing  array.  The  targeted  spherical  imaging  region  (18-cm  diameter)  is  depicted  at  isocenter,  d:  End-view  photo  of  the  Halbach  magnet 
mounted  on  high  friction  rollers.  Magnet  was  constructed  with  ABS  plastic  and  square  fiberglass  tubes  containing  the  NdFeB  magnets. 
Faraday  cage  not  shown. 


covered  with  a  high  friction  urethane.  The  MRI  console 
is  used  to  drive  a  stepper  motor  (model  34Y106S-LW8, 
Anaheim  Automation,  Anaheim,  CA,  USA)  that  is 
attached  to  the  aluminum  axle  of  one  of  the  rollers 
through  a  5:1  gearbox  (model  GBPH-0901-NS-005,  Ana¬ 
heim  Automation).  Magnet  rotation  is  incorporated  into 
the  pulse  sequence  so  that  it  is  controlled  by  the  MRI 
console  to  a  precision  of  one  degree  at  a  rate  of  10  deg/s. 
Peripheral  nerve  stimulation  is  not  a  concern  with  this 
B0  rotation  rate.  Even  at  10 x  the  current  rotation  rate, 
the  dB/dt  from  the  rotating  magnet  is  two  orders  of  mag¬ 
nitude  below  the  dB/dt  generated  by  a  modest  clinical 
gradient  system.  The  magnet  assembly  is  enclosed  in  a 
copper  mesh  Faraday  cage  to  reduce  RF  interference. 

An  initial  3D  field  map  was  obtained  with  a  three-axis 
gaussmeter  probe  attached  to  a  motorized  stage.  The 
measured  field  shape  is  roughly  quadrupolar,  similar  to 
the  fields  used  in  the  initial  realization  of  multipolar  Pat- 


Loc  (Parallel  Imaging  Technique  using  Localized  Gra¬ 
dients)  encoding  (13),  but  with  significant  higher-order 
components  as  well.  The  measured  field  variation  range 
in  y-z  (imaging  plane),  x-z,  and  x-y  planes  of  a  16-cm 
sphere  were  Afyz  =  95  kHz,  Afxz  =  60  kHz  and  Afxy  =  52 
kHz.  Large  Larmor  frequency  bandwidths  make  it  diffi¬ 
cult  to  design  RF  excitation  and  refocusing  pulses  that 
achieve  the  same  flip  angle  and  phase  across  all  the 
spins.  In  addition,  it  is  difficult  to  make  transmit  and 
receive  coils  uniformly  sensitive  over  the  entire  band¬ 
width.  Therefore,  shimming  was  done  to  decrease  field 
variation  (no  attempt  was  made  to  reshape  the  SEM). 
The  field  variation  was  shimmed  down  to  Afyz  =  32  kHz, 
Afxz  =  32.5  kHz  and  Afxy  =  19  kHz  with  the  addition  of 
small  shim  magnets  (0.5"  diameter,  0.25"  length  cylindri¬ 
cal  NdFeB  magnets)  which  were  attached  to  the  fiber¬ 
glass  rungs.  The  three  planes  of  the  shimmed  field  map 
are  shown  in  Figure  2. 
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FIG.  2.  Gaussmeter- measured  field  maps  of  the  spatial  encoding  magnetic  field  (SEM)  in  the  z-y  (imaging  plane),  z-x,  and  y-x  planes  of 
shimmed  the  Halbach  magnet.  The  B0  field  is  oriented  in  the  z  direction. 


An  accurate  field  map  is  critical  for  image  reconstruc¬ 
tion,  particularly  when  nonlinear  encoding  fields  force 
the  use  of  iterative  matrix  solvers  rather  than  the  Fourier 
transform  (14).  The  field  is  perturbed  by  external  fields 
(including  the  earth’s  magnetic  field),  and  must  be 
remapped  when  the  magnet  is  relocated.  To  quickly 
acquire  center-plane  field  maps,  a  linear  array  of  seven 
field  probes  spaced  1.5  cm  apart  was  constructed  (Fig. 
3a).  The  field  probes  are  tuned  5-mm-long,  4-mm-diame- 
ter,  18-turn  solenoids  measuring  signal  from  1-mm  capil¬ 
laries  of  CuS04-doped  water  (23).  To  acquire  a  field  map 
the  probes  are  held  stationary  while  the  magnet  is 
rotated  around  them.  Polynomial  basis  functions  are 
then  fit  to  the  measured  points  and  the  field  map  (Fig. 
3b)  is  synthesized.  The  polynomial  coefficients  up  to  6th 
order  of  one  magnet  rotation  angle  are  shown  in  Table  1. 
The  net  magnetic  field  from  the  NdFeB  magnets  is  sensi¬ 
tive  to  temperature  (on  the  order  of  4  kHz/deg  C  for  the 
Halbach  magnet)  as  well  as  interactions  with  external 
fields,  so  an  additional  field  probe  is  used  to  monitor 
field  drift  during  data  acquistion.  This  navigator  probe  is 
mounted  to  the  Halbach  array  and  rotates  with  the  mag¬ 
net.  The  measured  field  changes,  AB0,  are  then 
accounted  for  in  the  image  reconstruction. 

Acquisition  Method 

To  acquire  data,  the  magnet  is  physically  rotated  around 
the  sample  in  discrete  steps.  At  each  rotation  step,  gener¬ 
alized  projections  onto  the  nonlinear  field  are  acquired 
(similar  to  those  described  in  Schultz  et  al)  (24).  Exam¬ 


ples  of  these  projections  are  shown  in  Figure  4  for  a  sim¬ 
ple  two-sphere  phantom.  The  field  experienced  by  the 
spheres  changes  at  each  rotation  due  to  the  nonlinear 
SEM,  providing  new  information  in  each  projection. 

The  constructed  Rx  coil  array  (Fig.  5a)  consists  of  eight 
8-cm-diameter  loops  of  wire  encircling  the  FOV  on  the 
surface  of  a  14-cm-diameter  cylinder.  The  inductances  of 
the  coils  are  roughly  230  nH,  requiring  capacitors  on  the 
order  of  10  nF  (Voltronics,  Salisbury,  MD)  for  tuning. 
Geometric  decoupling  and  PIN  diode  detuning  were 
implemented  (25,26).  The  coils  are  tuned  and  matched 
to  50  Ohm  impedance  low  noise  preamplifiers  (MITEQ 
P/N  AU-1583,  Hauppauge,  NY,  USA). 

A  Tecmag  Apollo  console  with  TNMR  software  (Houston, 
TX,  USA)  was  used.  The  console  has  one  transmit  channel, 
three  gradient  channels,  and  one  receive  channel.  Because 
the  programmable  gradient  analog  outputs  are  not  needed  for 
gradient  coils,  they  are  used  for  other  purposes.  For  example, 
the  Gz  gradient  output  is  used  to  control  the  stepper  motor 
for  magnet  rotation.  The  fact  that  the  console  only  has  1 
receive  channel  means  that  true  parallel  imaging  cannot  be 
performed.  Instead,  the  receive  channel  is  switched  between 
the  coils  in  the  array,  acquiring  data  serially.  The  Gx  gradient 
output  along  with  a  RelComm  Technologies  (Salisbury,  MD, 
USA)  relay  and  Arduino  UNO  board  are  used  to  switch 
between  the  receive  coils.  Although  preamp  decoupling  has 
not  been  implemented  yet,  data  is  being  acquired  from  one 
coil  at  a  time,  permitting  the  other  receive  coils  to  be  detuned 
to  prevent  coupling. 

Two  scanner  coordinate  systems  are  defined  because 
the  object  and  RF  coils  remain  stationary  while  B0  is 


Table  1 

Polynomial  Composition  of  Halbach  Magnet  SEM  in  Hz/cm(m+n)a 
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z°y°  :  3.3e6 
zV  :  -62 
z°y2  :  164 
z°y3  :  3.9 
z°y4  :  0.95 
z°y5  :  -3.9e-3 
z°y6  :  9.2e-3 

zY  :  -89 
zV  :  104 
zY  :  -13.3 
zY  :  6.5 
zY  :  0.21 
zY  :  -6.6e-2 

z2y°  :  -274 
z2y1  :  -8.3 
z2y2  :  -0.53 
z2y3  :  4.4e-2 
z2y4  :  -1 ,9e-2 

z3y°  :  1 .9 
zY  :  -1-7 
z3y2  :  -0.12 
z3y3  :  -2.3e-2 

z4y°  :  1.1 e-2 
zY  :  -4.6e-2 
z4y2  :  0.11 

z®y°  :  2.4e-2 
zY  :  -1.8e-3 

z®y°  :  9.2e-3 

aThe  calculated  polynomial  coefficients  composing  the  z-y  plane  (2D  imaging  plane)  of  the  Halbach  spatial  encoding  field  are  shown. 
Measured  points  from  the  linear  array  of  field  probes  (Fig.  3)  were  used  for  this  sixth  order  polynomial  fit. 
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FIG.  3.  a:  Linear  array  of  seven  NMR  field  probes  used  for  mapping  the  static  magnetic  field.  The  probes  are  held  stationary,  while  the 
magnet  is  rotated  around  them  and  points  on  the  2D  center  plane  are  sampled,  b:  Measured  field  map  for  the  center  transverse  slice 
through  the  magnet  after  fitting  6th  order  polynomials  to  the  probe  data.  The  black  dots  mark  the  location  of  the  probe  measurements. 
The  field  is  plotted  in  MHz  (proton  Larmor  frequency).  This  field  distribution  serves  as  the  SEM  information  used  in  image 
reconstruction. 


rotated.  The  rotating  coordinate  system  of  the  magnet 
and  the  spins  is  defined  as  x7,  y7,  z7  (examples  shown  in 
Fig.  6b),  and  the  fixed  coordinate  system  for  the  coils 
and  objects  is  defined  as  x,  y,  z  (shown  in  Figs.  1,2). 
Image  reconstruction  requires  accurate  knowledge  of 
the  coil  sensitivity  map,Cgr(^)-  Here  the  index  q  refers 
to  the  coil  channel  and  r  to  the  rotation  position  of  the 
magnet.  The  coil  sensitivity  map  is  different  for  each 
rotation  position  because  is  formed  from  a  projec¬ 
tion  of  the  coil’s  Bt  field  onto  the  x'-y7  plane  (which 
rotates  with  the  magnet).  In  conventional  MRI,  B2~  is 
mapped  by  imaging  a  phantom  with  fully  sampled 
encoding  by  the  gradient  waveforms.  However,  this 
approach  is  not  possible  with  our  encoding  scheme 
because  knowledge  of  Cgr(x)  is  necessary  to  form  an 
image  without  aliasing. 

Because  of  the  difficulty  of  measuring  B on  our 
scanner,  we  use  estimated  B2~  maps.  Magnetostatic 
approximations  are  suitable  at  the  3.29  MHz  Larmor  fre¬ 
quency,  so  Bt  of  the  individual  coils  was  modeled  with 
Biot-Savart  calculations.  By  symmetry,  the  x  component 
of  the  circular  surface  coils’  B2  is  zero  in  the  center 
plane  FOV.  The  x7  component  of  B0  is  also  nearly  zero 
because  of  the  geometry  of  the  magnet.  So  the  coil  sensi¬ 
tivity  calculation  reduces  to  a  two  dimensional  problem, 


because  only  the  B2  component  perpendicular  to  B0  con¬ 
tributes  to  the  sensitivity  map. 

To  calculate  the  coil  sensitivity  map  for  each  rotation 
(r),  the  B2  component  parallel  to  B0r  (the  B0  vector  for 
rotation  r)  is  subtracted  and  we  are  left  with  the  perpen¬ 
dicular  component. 

B±t  =  B1-(B1-B0r)B0r.  [1] 

The  phase  is  equal  to  the  angle,  0r,  between  B^r  and  B0r, 
which  will  either  be  +90°  or  —90°  due  to  the  symmetry 
properties  discussed  above.  The  variation  in  a  single  coil’s 
B{~  as  a  function  of  B0  angle  is  illustrated  in  Figure  6a, b, 
and  the  Bt~  magnitude  for  four  coils  and  a  single  B0  angle 
is  shown  in  Figure  6c.  When  B0  points  along  the  normal 
to  the  coil,  the  sensitivity  profile  resembles  a  “donut”  pat¬ 
tern  with  low  sensitivity  in  the  center  of  the  FOV.  Maxi¬ 
mum  signal  sensitivity  occurs  when  B0  is  oriented 
orthogonal  to  the  normal  vector  of  the  coil  loop. 

Similar  to  single-sided  imaging  methods  (7),  echo  for¬ 
mation  requires  the  use  of  spin  echo  sequences  in  the 
presence  of  the  inhomogeneous  field.  The  T2*  of  the  sig¬ 
nal  is  short  due  to  the  static  SEM  and  it  is  impossible  to 
do  the  equivalent  of  gradient  echo  refocusing  because  the 
sign  of  the  SEMs  cannot  be  quickly  switched.  However, 


FIG.  4.  Schematic  depiction  of 
the  generalized  projections  (bot¬ 
tom  row)  of  an  object  onto  the 
rotating  SEM  field.  The  object 
consists  of  two  water-filled 
spheres  depicted  as  dashed 
black  lines  which  are  superim¬ 
posed  on  the  Halbach  magnet’s 
SEM  field  at  a  few  rotations 
(black  arrow  depicts  B0  orienta¬ 
tion).  The  NMR  spectrum  was 
acquired  with  a  single  volume  Rx 
coil. 
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FIG.  5.  a:  Photo  of  the  eight-channel  receiver  array  coil  with  3D  printed  disk-phantom  at  isocenter.  The  14-cm-diameter  array  is  made 
up  of  eight,  8-cm  loops  overlapped  to  reduce  mutual  inductance,  b:  Relative  voxel  size  is  illustrated  as  a  function  of  radius  from  the 
center  using  two  rotations  of  the  magnet’s  SEM  (field  isocontour  lines  illustrated  in  figure).  Symmetry  of  the  isocontours  causes  aliasing 
of  each  voxel  through  the  origin.  Using  the  local  sensitivity  profiles  of  an  encircling  array  of  coils,  the  correct  location  of  each  signal 


source  in  the  FOV  can  be  resolved,  (adapted  from  (28))  (C)  Photo  of 

the  encoding  can  be  repeated  and  averaged  to  improve 
SNR  in  a  spin-echo  train,  which  does  refocus  the  SEM. 
Unlike  in  high-field  systems,  the  specific  absorption  rate 
(SAR)  from  the  consecutive  180°  pulses  is  negligible 
because  of  the  low  excitation  frequency  (3.29  MHz). 

Unlike  conventional  MRI  scanners,  the  B0  field  of  the 
Halbach  magnet  is  oriented  radially  instead  of  along  the 


the  25  turn,  20-cm  diameter,  25-cm  length  solenoid  transmit  coil. 

bore  of  the  magnet.  This  means  that  in  order  for  B2+  to 
be  orthogonal  to  B0  at  all  rotations,  it  should  be  directed 
along  the  cylindrical  axis  of  the  Halbach  magnet.  This 
makes  a  solenoid  more  suitable  than  a  birdcage  coil  for 
radiofrequency  (RF)  excitation.  The  constructed  sole¬ 
noid,  shown  in  Figure  5c,  has  a  20-cm  diameter  and  a 
25-cm  length.  N  =  25  turns  of  AWG  20  was  chosen  as  a 


FIG.  6.  Biot-Savart  calculation  of  the  sensitivity  map  of  the  Rx  coil  array.  The  white  arrows  show  representative  orientations  of  B0,  which 
define  the  spin  coordinate  system  orientation  (x',  y',  z').  Image  reconstruction  requires  accurate  coil  sensitivity  profiles  for  each  B0  angle 
used  in  the  experiment.  a,b:  B-,~  magnitude  and  phase  for  a  single  representative  surface  coil  located  at  the  right  side  of  the  FOV  (position 
marked  with  white  line).  Because  of  the  symmetry  of  the  coils’  at  isocenter,  the  coils’  x'  component  is  approximately  zero,  and  the  process 
of  taking  the  projection  onto  thex'-y'  plane  (to  solve  for  B-,~)  will  produce  a  vector  parallel  or  antiparallel  to  y'.  Therefore,  the  B{~  phase  is 
always  +90°  or  -90°  in  the  depicted  transverse  isocenter  plane,  c:  Bf  magnitude  of  four  different  coils  of  the  array  (marked  with  white 
lines)  for  a  single  magnet  rotation  position. 
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reasonable  value  in  the  tradeoff  between  B2+  homogene¬ 
ity  and  parasitic  capacitance  from  closely  spaced  wind¬ 
ings.  The  70  uH  Tx  coil  is  tuned  to  3.29  MHz  with  eight 
230  pF  series  capacitors  distributed  along  the  length  of 
the  solenoid,  which  reduces  the  susceptibility  to  stray 
capacitance.  Because  the  static  SEM  field  is  “always  on”, 
the  transmit  coil  must  have  a  relatively  low  Q  to  excite  a 
wide  bandwidth  of  spins.  The  Q  of  the  coil  is  approxi¬ 
mately  60  corresponding  to  a  55  KHz  bandwidth.  A  1 
KW  power  amplifier  (Tomco,  Stepney,  SA,  Australia)  is 
used  to  produce  short  600  W  pulses  for  broadband  exci¬ 
tation  (25  |xs  for  90°  pulses  and  50  |xs  for  180°  pulses). 

PIN  diode  detuning  is  used  in  the  transmit  and 
receiver  coils  to  prevent  coil  interaction  (25).  The  tun¬ 
ing/matching  circuits  are  constructed  so  that  the  transmit 
coil  is  tuned  and  the  receive  coils  are  detuned  when  the 
pin  diodes  are  forward  biased  with  console  controlled 
DC  voltage.  The  converse  is  true  when  the  diodes  are 
reverse  biased  (Tx  coil  detuned  and  Rx  coils  tuned). 


Reconstruction  Method 

The  Halbach  magnet’s  spatial  encoding  field  is  approxi¬ 
mately  quadrupolar  and,  therefore,  produces  a  nonbijec- 
tive  mapping  between  object  space  and  encoding  space. 
This  encoding  ambiguity  leads  to  aliasing  in  the  image 
through  the  origin.  As  described  by  Hennig  et  al  (13), 
parallel  imaging  with  encircling  receive  coils  can  be 
used  to  disambiguate  the  nonbijective  mapping.  This  is 
possible  because  the  coil  sensitivity  profiles  provide 
additional  spatial  encoding  that  localizes  signal  within 
each  source  quadrant  of  the  FOV,  eliminating  aliasing. 
This  idea  is  illustrated  in  Figure  5b.  This  specific  imple¬ 
mentation  of  the  portable  scanner  closely  resembles  the 
case  of  PatLoc  imaging  with  quadrupolar  fields  and  a 
radial  frequency-domain  trajectory  (28).  However,  the 
measured  Halbach  SEM  is  not  purely  quadrupolar,  and 
the  presence  of  arbitrary  field  components  prevents  the 
decomposition  of  the  rotating  SEM  into  linear  combina¬ 
tions  of  two  orthogonal  encoding  fields.  For  this  reason, 
the  direct  back-projection  reconstruction  method 
described  in  Schultz  et  al  (28)  is  not  valid,  and  iterative 
matrix  methods  such  as  those  described  in  Stockmann 
et  al  (29)  are  used. 

The  discretized  signal  acquired  by  a  coil  [q]  at  a  given 
magnet  rotation  (r)  at  time  n  can  be  described  as 

sq,r(n)  =  J2  Cq,r(x)e-i2l,k{r*’n)m(x)  [2] 

X 

where  m(x)  is  the  magnetization  of  the  object  at  location 
vector  x,  Cgr(x)  is  the  complex  sensitivity  of  the  coil  at 
location  x,  and  Jc(r,  x,  n)  is  the  evolved  phase  from  the  non¬ 
linear  gradient  at  rotation  r,  location  x,  and  time  n.  The 
exponential  term  and  coil  sensitivity  term  can  be  grouped 
together  to  form  the  encoding  function  encq:r(x ,  n ). 

Sq,r(n)  =  'Ll  encv{x,  n)m(x)  [3] 

X 

The  matrix  form  of  this  signal  equation  for  a  single 
projection  readout  acquired  with  one  RF  coil  is 
simply 


Sqr  —  Egrim  [4] 

The  acquired  signal,  Sq>r,  is  a  vector  made  up  of  the 
sampled  readout  points  (Nsmp).  The  object  that  we  are 
solving  for,  m,  is  a  vector  made  up  of  all  the  image  vox¬ 
els  (Nvox).  The  encoding  matrix,  Eqr,  contains  the 
evolved  phase  of  each  voxel  in  the  FOV  for  each  time 
point  in  the  acquisition  as  well  as  the  coil  sensitivity 
multiplier.  With  linear  gradient  fields,  E  is  made  up  of 
the  sinusoidal  Fourier  basis  set,  which  allows  the  image 
to  be  reconstructed  using  radial  back-projection,  k-space 
re-gridding,  and  other  approaches.  In  the  nonlinear  SEM 
case,  E  is  more  complicated,  but  can  be  calculated  from 
the  measured  field  maps.  Before  the  appropriately 
rotated  field  map  is  used  to  calculate  the  phase  evolu¬ 
tion,  the  field  change  captured  by  the  navigator  probe 
during  the  acquisition  is  added  as  a  global  offset. 

A  separate  block  of  the  encoding  matrix,  Eq?r  is  calcu¬ 
lated  for  the  data  acquired  by  each  coil  at  each  rotation. 
There  will  be  a  total  of  R*C  blocks  (where  R  is  total 
number  of  rotation  and  C  is  the  total  number  of  coils), 
which  are  vertically  concatenated  to  form  the  full  encod¬ 
ing  matrix,  E.  S  is  also  made  up  of  vertically  concaten¬ 
ated  subparts,  Sqr,  which  are  the  signals  acquired  from 
each  coil  at  each  rotation. 

To  reconstruct  the  image  from  the  acquired  data,  the 
object,  m,  can  be  found  by  inverting  the  matrix,  E.  How¬ 
ever,  the  full  encoding  matrix  size  is  Nsmp  *R*C  x  Nvox. 
In  the  typical  case  of  256  readout  points,  181  rotations, 
eight  coils,  and  a  256  x  256  voxel  reconstructed  image, 
the  full  matrix  size  is  371K  x  65K.  Because  it  is  not  com¬ 
putationally  feasible  to  invert  this  matrix,  iterative  meth¬ 
ods  such  as  the  Conjugate  Gradient  method  (30)  and  the 
Algebraic  Reconstruction  Technique  (31,32)  can  be  used 
to  solve  for  the  minimum  norm  least  squares  estimator  of 
m.  The  generality  of  this  approach  allows  arbitrary  field 
shapes  and  coil  profiles  as  well  as  systematic  errors  such 
as  temperature-dependent  field  drifts  to  be  incorporated 
into  the  encoding  matrix. 

The  reconstructed  images  and  simulations  shown  here 
were  done  using  the  Algebraic  Reconstruction  Tech¬ 
nique.  The  encoding  matrix  was  calculated  line  by  line 
during  the  reconstruction  using  the  appropriately  rotated 
and  temperature  drift-corrected  field  map  and  the  calcu¬ 
lated  coil  sensitivity  profiles  for  the  given  B0  direction. 
To  demonstrate  the  importance  of  temperature  drift 
compensation,  a  phantom  image  was  also  reconstructed 
with  an  uncorrected  encoding  matrix.  The  field  of  view 
of  the  images  is  16  cm  and  the  in-plane  voxel  size  is 
0.625  mm. 

Phantom  Imaging  Methods 

Images  of  a  “MIT/MGH”  phantom  were  acquired  both 
with  a  single  channel  solenoid  Rx  coil  and  with  seven 
coils  of  the  Rx  array.  The  3D  printed  polycarbonate 
phantom  is  1.7-cm  thick  with  a  13-cm  diameter,  and  is 
filled  with  CuS04-doped  water.  Thirty-two  averages  of  a 
six  spin-echo  train  (TR  =  550  ms,  echo-spacing  =  8ms) 
were  acquired  for  91  rotation  angles  over  180  degrees. 
Navigator  field  probe  data  was  also  acquired  at  each  rota¬ 
tion.  The  coil  array’s  lengthy  acquisition  time  of  66  min 


43 


8 


Cooley  et  al. 


f#  t  ^ 

la  ^ 

ll 

a 

b 

c 

FIG.  7.  Experimental  256  x  256  voxel,  16-cm  FOV  images  of  a  3D  printed  phantom  with  CuS04  doped  water  occupying  the  interior  of 
the  letters  and  polycarbonate  plastic  surrounding  it.  The  phantom  has  a  13-cm  diameter  and  is  1.7  cm  thick.  91  magnet  rotations 
spaced  2°  apart  were  used,  readout  bandwidth/Npts  =  40  KHz/256,  TR  =  550  ms,  spin-echo  train  length  =  6  or  16,  with  8  ms  echo¬ 
spacing.  Echoes  in  the  spin-echo  train  for  a  given  rotation  were  averaged,  a:  Image  acquired  with  solenoid  Rx  coil  (32  averages  of  a  6 
spin-echo  train),  b:  Image  acquired  with  seven  coils  of  the  Rx  coil  array  (8  averages  of  a  16  spin-echo  train).  Temperature  drift  was  not 
corrected  for.  c:  Image  from  same  data  as  (B),  but  with  temperature  drift  correction  implemented. 


results  from  multiplexing  a  single  console  receiver  and 
would  be  reduced  to  7.3  min  by  acquiring  data  from  all 
channels  and  the  field  probe  in  parallel. 

A  1-cm-thick  lemon  slice  was  imaged  using  only  the 
bottom  five  surface  coils  with  181  1°  rotations.  The  total 
acquisition  time  was  93  minutes  (15.5  min  if  surface 
coils  and  navigator  probe  were  acquired  in  parallel).  A 
single  average  of  a  128  echo  train  at  each  rotation  pro¬ 
vided  sufficient  SNR.  Each  echo  was  recorded  as  256  pts 
with  a  40  KHz  BW  (TR  =  4500  ms,  echo  spacing  =  8 
ms).  For  comparison,  the  lemon  image  was  also  recon¬ 
structed  using  only  91  rotations  of  181  acquired  rotations 
in  addition  to  the  full  reconstruction. 

Image  Simulation  Methods 

The  described  acquisition  method  was  simulated  using 
the  measured  field  map  from  the  central  slice  of  the  Hal- 
bach  magnet.  Images  were  simulated  using  a  high  resolu¬ 
tion  Ta  brain  image  or  a  numerically  generated 
checkerboard  with  2.5-mm  grid  size  as  the  “object.”  The 
measured  field  map  and  calculated  coil  profiles  of  the 
eight  coil  array  were  used  in  the  forward  model  to  gener¬ 
ate  the  simulated  data.  In  one  simulation,  an  artificial 
field  map  was  used  to  simulate  the  addition  of  a  linear 
field  component  to  measured  SEM.  Complex  noise  was 
added  to  the  simulations  to  match  noise  levels  observed 
in  comparable  phantom  projection.  These  simulations 
were  done  with  the  same  sequence  parameters  of  the 
lemon  image:  181  1°  magnet  rotations,  256  pt  readout,  40 
KHz  BW,  echo  spacing  =  8  ms. 

RESULTS 

Experimental  images  of  the  3D  printed  “MIT/MGH” 
phantom  are  shown  in  Figure  7.  The  image  acquired 
with  the  solenoid  coil  used  in  transmit/receive  mode  is 
shown  in  Figure  7a.  Only  the  “MGH”  part  of  the  phan¬ 
tom  was  filled  at  the  time,  so  the  top  half  of  the  image 
should  ideally  be  empty.  Instead  the  expected  aliasing 
pattern  is  seen  through  the  center  onto  frequency 


matched  quadrants  of  the  FOV.  The  aliased  image  is 
markedly  more  blurry  than  one  would  expect  for  a 
purely  quadrupolar  field,  which  maps  all  signals  sym¬ 
metrically  about  the  center  during  reconstruction.  This 
discrepancy  likely  arises  due  to  the  presence  of  first- 
order  and  higher-order  field  terms  which  perturb  the 
symmetry  of  the  dominant  quadrupolar  field. 

The  importance  of  monitoring  and  correcting  for  field 
drift  due  to  temperature  is  emphasized  by  comparing 
Figure  7b  and  Figure  7c  which  show  images  with  and 
without  temperature  drift  correction.  The  drift  correction 
is  achieved  by  monitoring  the  frequency  seen  by  the  nav¬ 
igator  field  probe  which  rotates  with  the  magnet.  This 
probe’s  frequency  is  ideally  independent  of  the  rotation 
angle  during  the  acquisition,  but  varies  due  to  two 
causes.  First,  small  changes  in  room  temperature  trans¬ 
late  to  a  global  scaling  of  the  Halbach  array’s  magnetiza¬ 
tion  and  thus  the  central  B0.  With  no  attempt  to  insulate 
or  stabilize  the  magnet’s  temperature,  changes  up  to  0.4 
°C  and  1.6  KHz  were  observed  over  an  hour.  The  second 
cause  for  the  fixed  probe’s  change  in  field  as  a  function 
of  rotation  is  due  to  the  changing  vector  sum  of  the 
earth’s  field  and  the  Halbach  field.  This  effect  creates  a 
peak-peak  variation  of  3.7  kHz  for  the  magnet  location 
and  orientation.  This  effect  must  also  be  incorporated  in 
the  encoding  matrix.  Even  though  field  drift  correction  is 
applied  to  Figure  7c,  some  of  the  letters  are  sharper  than 
others;  this  is  likely  attributable  to  field  map 
inaccuracies. 

The  1-cm  lemon  slice  images  are  seen  in  Figure  8.  The 
use  of  five  of  eight  coils  of  the  receive  array  prevents  ali¬ 
asing  in  the  image,  but  center  blurring  is  more  pro¬ 
nounced  in  these  images  than  in  the  simulations  (Fig.  9). 
Figure  8a  was  reconstructed  using  half  of  the  rotations 
angles  of  Figure  8b,  resulting  in  poorer  image  quality 
and  streaking  artifacts. 

Figure  9  shows  an  encoding  and  reconstruction  simu¬ 
lation  using  a  typical  high  field  Ta-weighted  brain  MRI 
as  the  imaging  object  (Fig.  9a).  Noise  was  added  to  the 
object  model  to  simulate  the  lower  SNR  of  the  low  field 
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FIG.  8.  Experimental  256  x  256 
voxel,  16-cm  FOV  image  of  a  1- 
cm  thick  slice  of  lemon  placed 
off  axis  in  the  magnet.  Five 
receiver  coils  of  the  array  were 
used  to  acquire  one  average  of 
a  128  spin-echo  train,  readout 
bandwidth/Npts  =  40  KHz/256, 
TR  =  4500  ms,  echo-spacing  = 
8  ms.  a:  91  magnet  rotations 
spaced  2°  apart  were  used  (b) 
181  magnet  rotations  spaced  1° 
apart  were  used. 


scanner.  Figure  9b  shows  a  simulated  image  using  the 
measured  encoding  field  of  Halbach  magnet.  There  is  no 
aliasing  in  the  image  because  the  calculated  coil  sensitiv¬ 
ity  profiles  of  the  eight  channel  Rx  array  were  used. 
However,  there  is  blurring  in  the  center  which  coincides 
with  the  shallow  region  of  the  nonlinear  gradient  field. 
The  center  blurring  in  Figure  9c  is  reduced  because  the 
image  was  simulated  using  an  artificial  field  map  that 
consists  of  our  measured  SEM  plus  an  additional  linear 
field  of  500  Hz/cm.  The  simulation  of  the  2.5-mm  grid 
numerical  phantom  (Fig.  9d)  shows  the  ultimate  resolu¬ 
tion  possible  with  the  existing  experimental  protocol  in 
the  absence  of  systematic  errors.  Outstanding  resolution 
at  the  periphery  gradually  gives  way  to  a  blurry  central 
region. 

DISCUSSION 

As  expected,  the  nonbijective  mapping  of  the  Halbach 
magnet’s  SEM  results  in  aliasing.  Fortunately,  as 


described  in  Hennig  et  al  (13)  the  aliasing  is  resolved 
by  the  addition  of  a  multichannel  receive  array  with 
differing  spatial  profiles  and  an  appropriate  geometry. 
Because  the  Halbach  encoding  is  dominated  by  the 
quadrupolar  “PatLoc”  SEM,  the  system’s  spatially  vary¬ 
ing  voxel  size  changes  approximately  as  c/p  within  the 
FOV,  where  p  is  the  radius  and  the  constant  c  depends 
on  the  strength  of  the  SEM  and  the  length  of  the  read¬ 
out  (27).  This  means  that  our  Halbach  magnet  encoding 
field  results  in  higher  resolution  at  the  periphery  due 
to  the  uniform  nature  of  the  SEM  near  the  center  of 
the  FOV.  This  center  blurring  is  seen  in  both  the 
experimental  images  in  Figure  8  and  the  simulations  in 
Figure  9. 

While  we  did  not  attempt  to  control  the  precise  spheri¬ 
cal  harmonic  distribution  in  the  magnet  design,  future 
work  will  likely  benefit  from  shimming  the  magnet  to 
obtain  a  more  desirable  SEM.  For  example,  if  a  sufficient 
linear  term  were  added,  the  uniform  encoding  field  region 
would  not  lie  on-axis  with  the  rotation.  In  this  case, 


FIG.  9.  Simulated  images  using  the  calculated  sensitivity  profiles  of  the  eight  coil  Rx  array  to  generate  the  forward  model  for  181  1° 
rotations  of  the  encoding  field,  6.4  ms,  256  point  readouts.  The  data  seen  by  the  Halbach  scanner  was  simulated  by  processing  this 
“object”  through  the  forward  model  and  adding  noise  to  make  it  consistent  with  the  SNR  of  the  time-domain  signals  measured  in  a 
water  phantom.  The  model  data  was  then  reconstructed  using  the  Algebraic  Reconstruction  Technique  in  a  16-cm  FOV.  a:  Reference 
high  resolution  3T  T-i  weighted  brain  image  used  as  the  model  object.  Note:  the  SEMs  were  scaled  to  the  brain  FOV.  b:  Simulated 
reconstruction  using  the  measured  SEM  to  generate  the  forward  model,  c:  Simulated  reconstruction  using  the  measured  SEM  with  the 
additional  artificial  linear  field  component  (500  Hz/cm),  d:  Simulated  reconstruction  of  a  2.5-mm  grid  numerical  phantom.  Only  one 
quadrant  of  the  FOV  is  shown,  the  center  of  the  FOV  is  marked  with  white  cross-hairs  in  the  upper  right. 
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which  is  simulated  in  Figure  9c,  the  “blind-spot”  would 
move  around  the  object  allowing  some  rotations  to  con¬ 
tribute  to  encoding  of  any  given  pixel,  as  previously 
explored  in  “O-space  imaging”  (14).  Pursuing  this  strategy 
even  further  would  result  in  a  SEM  containing  only  a  lin¬ 
ear  term.  In  this  case,  the  encoding  becomes  very  similar 
to  a  radial  imaging  scheme  with  conventional  gradients, 
and  to  the  strategy  proposed  by  Cho  et  al.  who  used  a 
rotating  gradient  coil  in  a  conventional  magnet  (12).  With 
accurate  field  mapping  instrumentation  and  shimming 
software,  we  suspect  that  the  magnet  could  be  shimmed 
to  a  more  desirable  SEM.  Although  a  linear  SEM  would 
eliminate  the  encoding  hole  and  allow  a  more  straightfor¬ 
ward  reconstruction  method,  there  are  advantages  to 
second-order  SEMs,  including  the  coincidence  of  the  high 
spatial  resolution  area  and  high  coil  sensitivity  area  near 
the  edge  of  the  FOV. 

The  lemon  images  of  Figure  8  show  that  when  91  pro¬ 
jection  rotations  are  used  instead  of  181,  a  radial  streak¬ 
ing  artifact  is  visible.  The  streaking  artifacts  are 
consistent  with  those  arising  in  conventional  under¬ 
sampled  radial  trajectories  played  by  linear  SEMs  as 
well  as  undersampled  radial  trajectories  played  by  Pat- 
Loc  SEMs  (27).  It  has  been  shown  that  the  use  of  total 
variation  and  total  generalized  variation  priors  during 
reconstruction  suppresses  streaking  artifacts  in  under¬ 
sampled  conventional  radial  (33)  and  PatLoc  radial  (34) 
acquisitions.  Similar  techniques  may  be  pursued  in 
future  work  to  suppress  streaking  in  images  obtained 
with  fewer  projection  rotations  of  our  scanner. 

The  simulations  in  Figure  9  show  the  theoretical 
resolution  of  the  scanner  when  systematic  errors  are 
eliminated.  These  errors  are  most  likely  a  result  of  field 
map  and  coil  sensitivity  profile  inaccuracies,  which  are 
critical  to  the  iterative  reconstruction  (14).  The  current 
coil  sensitivity  profiles  facilitate  proof-of-concept 
reconstructions,  but  their  fidelity  is  suspect  because 
they  were  calculated  rather  than  measured.  In  these  cal¬ 
culations  the  magnetostatic  Biot  Savart  approximation 
was  used  with  no  external  structures  present.  While 
wavelength  effects  in  the  body  are  not  expected  at  this 
frequency,  the  close  proximity  of  the  conducting  mag¬ 
nets  and  other  coils  might  perturb  the  experimental 
fields.  Additionally,  a  2D  field  map  is  currently  used  to 
reconstruct  thin  samples  (1  to  1.5  cm  thick),  but  field 
variation  does  exist  in  the  x  direction  (along  the  axis  of 
the  Halbach  cylinder)  within  the  sample  thickness. 
This  causes  through-plane  dephasing  and  must  be 
incorporated  into  the  encoding  matrix  based  on  a  3D 
field  map. 

Field  map  errors  arise  from  temperature  drifts  which 
are  significant  on  the  time  scale  of  the  imaging  and  map¬ 
ping  acquisitions.  We  have  shown  that  any  uncorrected 
temperature  drift  causes  substantial  blurring  in  the  image 
(Fig.  7b).  Temperature  drift  is  a  pervasive  problem  in 
permanent  magnet  MRI  and  has  been  addressed  in  sev¬ 
eral  ways.  In  the  current  experimental  protocol  the  fre¬ 
quency  at  a  fixed  point  is  measured  at  every  rotation  and 
the  drift  is  built  into  the  encoding  matrix  as  a  global  off¬ 
set  to  the  field  maps.  This  method  reduces  blurring  con¬ 
siderably  (Fig.  7c),  but  other  options  have  been  proposed 
for  permanent  magnet  NMR  and  MRI  that  may  offer 


higher  encoding  matrix  accuracy.  For  example,  Rose  and 
Haishi  (1)  describe  the  implementation  of  a  NMR  lock 
method  plus  thermal  insulation.  Additionally,  a  new 
Halbach  design  was  recently  reported  which  uses  two 
types  of  magnet  materials  with  different  temperature 
coefficients  to  substantially  reduce  the  effect  of  tempera¬ 
ture  changes  on  the  B0  field  (35).  When  compared  with 
an  uncompensated  SmCo  magnet,  this  reduced  their  tem¬ 
perature  coefficient  100-fold,  bringing  the  field  drift 
down  to  10  ppm  for  a  3°C  temperature  change.  However 
this  method  has  the  disadvantage  of  producing  a  lower 
field  and  requiring  more  magnet  material  than  the  tradi¬ 
tional  design. 

For  time-efficient  acquisitions,  true  parallel  imaging 
will  be  needed.  To  accomplish  this  goal,  a  multichannel 
receiver  console  is  required,  as  well  as  the  implementa¬ 
tion  of  preamp  decoupling.  This  is  advantageous  for 
practical  diagnostic  reasons  and  will  also  alleviate  the 
field  drift  problem  by  shortening  acquisition  times.  In 
addition  to  multiple  channels,  future  prototypes  must  be 
made  larger  to  accommodate  the  human  head.  Although 
the  head  can  be  fit  into  the  presented  magnet,  its  36-cm 
diameter  does  not  leave  sufficient  room  for  the  transmit 
and  receive  arrays  as  well  as  the  structural  supports  for 
the  magnetic  material.  Construction  of  a  larger  diameter 
magnet  with  the  same  basic  design  will  result  in  a 
reduced  B0  field,  although  this  could  be  mitigated  by 
adding  more  magnet  material  and/or  higher  grade  mate¬ 
rial.  The  current  B0  field  of  77.3  mT  is  estimated  to 
decrease  to  62  mT  if  the  diameter  is  increased  to  40  cm. 
However  if  24  N45  NdFeB  magnet  rungs  are  used  instead 
of  20  N42  rungs,  a  field  of  80  mT  is  theoretically  achiev¬ 
able.  The  standard  landmark  for  brain  imaging  (between 
the  eyebrows)  is  18  cm  above  the  shoulders.  The  pre¬ 
sented  Halbach  magnet  was  designed  using  the  maxi¬ 
mum  cylinder  length  that  allows  the  brain  to  be  centered 
in  the  magnet  (2  x  18  cm).  Future  magnet  designs  will 
likely  adhere  to  this  constraint  because  increasing  the 
length  of  the  magnet  requires  increasing  the  bore  size  to 
fit  shoulders,  which  would  result  in  a  considerably 
weaker  B0. 

In  the  described  experiments  the  B0  field  rotates  rela¬ 
tive  to  the  receiver  coils  (coils  are  stationary),  which 
causes  the  shape  of  the  coil  profiles  to  change  with  each 
acquisition  angle.  However  this  arrangement  is  not  a 
requirement  for  rotating  SEM  imaging,  and  in  theory  the 
receiver  coils  could  rotate  with  the  magnet.  In  this  case, 
the  coil  sensitivity  profiles  are  simply  rotated  for  each 
acquisition  angle,  but  the  shapes  of  profiles  do  not 
change.  Data  acquisition  with  rotating  coils  and  station¬ 
ary  coils  was  simulated.  However,  there  was  not  a  signif¬ 
icant  difference  in  performance  in  either  the  visual 
appearance  of  the  reconstructed  images  or  the  RMSE 
(root  mean  squared  error).  For  data  simulated  with  91 
magnet  rotations  there  was  a  0.2%  RMSE  improvement 
when  using  the  rotating  coil  profiles,  and  for  data  simu¬ 
lated  with  23  magnet  rotations  (undersampled)  there  was 
a  3.6%  RMSE  improvement.  This  suggests  that  rotating 
the  coil  array  may  improve  performance  when  data  is 
undersampled.  The  23  rotation  simulated  images  are 
included  in  the  supplemental  material.  The  rotating 
receive  coil  case  is  similar  to  the  RRFC  (Rotating  RF 
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Coils)  method  described  in  Trakic  et  al  and  Li  et  al 
(36,37),  where  continuously  rotating  surface  coils  are 
used  in  a  conventional  magnet  for  parallel  imaging. 

The  goal  of  the  current  work  was  to  provide  a  proof-of- 
principle  that  the  basic  2D  encoding  scheme  can  be  per¬ 
formed,  which  was  demonstrated  with  2D  imaging  of 
thin  samples.  However,  the  addition  of  third  axis  encod¬ 
ing  is  an  obvious  requirement  for  medical  applications. 
A  promising  possibility  for  encoding  the  third  dimension 
(along  the  axis  of  rotation)  is  TRASE  (15,16).  TRASE 
uses  custom-designed  RF  coils  to  generate  uniform 
amplitude  but  linear  B2+  phase  variation  along  the 
encoding  axis.  Spatial  encoding  is  achieved  using  at 
least  two  Tx  coils  with  different  phase  gradients  (typi¬ 
cally  differing  by  their  sign).  Spin-echo  trains  are  used 
in  which  the  linear  phase  variation  is  changed  by  180 
degrees  in  between  successive  refocusing  pulses.  As  the 
sign  of  the  refocusing  pulse  is  flipped  over  the  course  of 
the  echo  train,  k-space  is  traversed  one  echo  at  a  time. 
The  resolution  depends  on  the  number  of  echoes  used 
and  the  slope  of  the  transmitted  B2  phase  ramp  across 
the  FOV  (16).  The  approach  is  synergistic  with  the  echo 
trains  used  in  the  presented  encoding  scheme  for  pur¬ 
poses  of  signal  averaging.  Furthermore,  at  low  field, 
TRASE  spin-echo  trains  do  not  suffer  from  the  SAR  lim¬ 
its  that  may  impact  the  method’s  performance  at  high 
field. 

CONCLUSIONS 

Using  an  inhomogeneous  magnet  for  spatial  encoding  in 
lieu  of  gradient  coils,  we  have  constructed  and  demon¬ 
strated  a  lightweight  scanner  for  2D  MR  imaging  with 
minimal  power  requirements.  The  2D  proof-of-concept 
images  from  this  nearly  head-sized  imager  show  the 
ability  of  this  encoding  scheme  to  produce  sufficient 
spatial  resolution  and  sensitivity  for  the  detection  and 
characterization  of  many  common  neurological  disorders 
such  as  hydrocephalus  and  traumatic  space-occupying 
hemorrhages.  Future  work  in  perfecting  the  calibration 
methods  is  likely  to  bring  experimental  image  quality 
closer  to  the  theoretical  limit,  but  the  resolution  of  the 
current  system  is  sufficient  for  identifying  gross  patholo¬ 
gies.  With  the  future  implementation  of  true  parallel 
imaging  and  3D  encoding,  this  scanner  has  the  potential 
to  enable  a  truly  portable,  low-cost  brain  imaging 
device. 
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